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(Aus dem Max Planck-Institut für Strömungsforsch ung, Göttingen) 


Für eine von E. Rei ßner angegebene und durchgeführte Verfeinerung der klassischen Plattenbiegungs- 

theorie werden die Grundgleichungen hergeleitet, bei-denen sich im Unterschied zur klassischen Theorie 

—- drei Randbedingungen befriedigen lassen. Als Anwendungsbeispiel wird die am Rand kin gelagerte Recht- 
eckplatte mit veränderlicher Belastung behandelt. 


Following the ideas of EB. Reissner, who refined the classical theory of the bending of plates, 


the author deduces the fundamental equations. Contrary to tue classical theory, Ihe method allows to salisfy 


re : ne conditions. In an Brompl, the rectangular plate with variable load and free boundaries 
is treate 


N RES Pour une amelioration de la theorie classique du courbement de plaques, laquelle fut donnee et r£alisee 

i par E. Reißner, les &qualions fondamentales sont deduites, au moyen desquelles, & la difference de la th£orie 
classique, trois conditions marginales peuvent Eire accomplies. Comme exemple est trailde la plaque 
rectangulaire avec des bords libres et avec une charge variable. 


Ana yıyymenna KiraccHgeckof TeOpuM IIPOTHÖAHHA TLIACTUH, MAHHOTO H IIPOBe/IEHHOTO 

9. Pelichepon, BEIBONATCH PyHNaMeHTanbHBIe yPaBHeHHg, Aa KOTOPEX, B IIPOTHBOTIONOHOCTE 

_ K KJIaCCHyecKOH TeOPHH, MO)KEeT OBITb YAOBIETBOpPeHO TpeM KpaeBbIM YcoBHfaMm. B kauectse 

UPHEIANIHOTO IIPmMepa PACcMaTpHBaeTcH CBODOAHAA IIO KpasnMm IIPAMOYTONBHAA ILIACTUHA_ 0 

tepemermoh HaTPy3KoH. 
AR 1. Einführung 

Der klassischen Biegungstheorie elastischer Platten, deren abschließende Darstellung 

| E Kirchhoff in seiner bedeutenden Abhandlung [1] vom Jahre 1850 gegeben hat, liegen be- 

_ kanntlich die beiden Voraussetzungen zu Grunde, daß einmal die Abmessung der Platte senk- 

_ recht zu ihrer Oberfläche, ihre Dicke, klein sei gegenüber den linearen Oberflächenabmessungen 

und zum anderen auch noch die Durchbie gung der Platte gering im Verhältnis zu ihrer Dicke. 

- Nun sind Plattenaufgaben nach der Kirchhoffschen Theorie zwar einer mathematischen 

Behandlung weitgehend zugänglich, jedoch kann dabei als Folge der ‘Vereinfachungen eine 

ernste grundsätzliche Schwierigkeit deswegen auftreten, weil sich nur zwei unabhängige Rand- 

bedingungen erfüllen lassen, während man ja im allgemeinen (z.B. bei der frei gelagerten Platte 

oder am freien Plattenrand) deren drei hat. Zur Behebung dieser Schwierigkeit hat Eric 


 Reißner vor einigen Jahren eine Verfeinerung ‚der klassischen Theorie angegeben, die ohne 


R Be uhermäßiee Steigerung des Rechenaufwandes zu einem erweiterten Differentiälgleichungssystem 
7 IE führt, bei dem es möglich und notwendig ist, sämtliche drei Randbedingungen zu befriedigen. 
Diese neue Theorie ist zweifellos schon in grundsätzlicher Hinsicht von Bedeutung, darüber 
hinaus aber auch von praktischem Nutzen, indem sie gewisse Fälle einzubeziehen gestattet, bei 

Bi funden stehen). 
N. - Da die grundlegenden A Drerächen Arbeiten [5], [6], [7], die kurz vor und nach Kriegs- 
ende erschienen sind, in Deutschland wohl bisher nicht allgemein bekannt geworden sind, so 
aR - dürfte eine kurze Darstellung der neuen Theorie für einen Teil des Leserkreises dieser Zeitschrift 


Re Theorie klar herauszuarbeiten, und müssen deshalb unmittelbar an die Ausgangsgleichungen der 

Elastostatik anknüpfen. Aber wir wollen uns nicht mit der Herleitung der erweiterten Grund- 

ee =  gleichungen begnügen, sondern auch ein Anwendungsbeispiel bringen, bei dem.die Reißnersche 

= Theorie in weitem Umfang zum Einsatz kommt, nämlich das Beispiel der Biegung einer am 
; es 2A A frei gelagerten Rechteckplatte unter veränderlicher Belastung. 


Sn 2% Grundlagen der verfeinerten Theorie 


au dem Innern der Platte (Dicke 2%) denke man sich ein qua 
das an 1 er ‚dem: N ist auch die res des zu- 


TR 


denen die Ergebnisse der klassischen Theorie nicht im Einklang mit den experimentellen Be- R 


- nicht ohne Interesse sein. “Wir sind dabei bestrebt, den Unterschied gegenüber der klassischen 


derförmiges Element 
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grunde gelegten cartesischen ee 2 Se ferner die gewohnte 
i sweise für die Spannungen in den Schnittebenen. R ; 
Se ringe Plattendicke legt es nahe, „über den Se = 
z-Richtung einfache Näherungsannahmen zu machen. Die Reißnersche Theorie . a 
diese aus der klassischen Theorie geläufigen Annahmen, die etwas später auch noch gen N 
muliert werden. Während nun aber dort die vertikalen Spannungskomponenten Oz, Tgz n Ei 
(in Bild 1 sind die beiden letzten durch die Strichstärke hervorgehoben) vernachlässigt (oder 


Bild 1 


höchstens nachträglich in die Rechnungeingeführt) werden, berücksichtigt die verfeinerte Theorie 
von vornherein auch diese Spannungskomponenten als gleichberechtigt neben den übrigen. 

In bekannter Weise werden aus den ‚Spannungen durch Zusammenfassung in Richtung 
der Dickenerstreckung auf die Längeneinheit bezogene ‚‚Schnittkräfte“: 


h h Be 
Q, — ta. d bzw. 9, = [tu.da nenne lab) 


gebildet sowie auf die Längeneinheit bezogene ‚‚Schnittmomente‘; und zwar liefern die Normal- 
spannungen die Biegemomente 


h N 
M,= [0o,: 2dz bzw. M= IK? Sg, ee (2a,b), 
op = 


und die horizontalen Schubspannungen das Drillmoment 


Bild 2 zeigt die Schnittkräfte und Schnittmomente für unser 
Plattenelement, ferner bedeutet p(z, y) die Belastung pro 


Flächeneinheit. 
Die Gleichgewichtsbedingungen 
tn =0 re (AR), 
ara —Q=0 . (4b), 
a & A (4c) 


ergeben zunächst drei Gleichungen für die fünf unbekannten 
Größen Q,, Q,; M,, M,, H. Die noch fehlenden Gleichun- 


Bild 2 gen liefern die Beziehungen, die zwischen den Verzer- _ 


rungen und Spannungen im Körper bestehen. 
Hierfür hat man nach dem verallgemeinerten Hookeschen Gesetz bei Annahme eines 
homogenen und isotropen Plattenmaterials?) und wegen der vorausgesetzten kleinen De- 
formationen die Gleichungen U ’ 


1 
&, = Fe = 51 —u(0y +0,)} ww OR ' . ee . .. (5a), 
?) Jedoch sei ausdrücklich darauf hingewiesen, daß die Reißnersche Theorie auch allgemeinere Fälle 


von nicht homogenem und nicht isotropem Plattenmaterial (wo also mehr als zwei unabhängige 
Elastizitätskonstanten auftreten) umfaßt (vgl. [7)). - E 
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ya u 0, + 0,)} Re en io), 
(== 40 — u (0, +0,)}) ee 
wege, „+0,)) RE. 
ou oV 1 
Pay 0Y * 0. Ge (5d), 
ou .9W 1 
ta2 7 7, ze . GQG": (8e), 
oV _0W 1 
Ywe=, ir 3 er N (Sf), 
E 
"rat 


Darin bedeuten e und y die Verzerrungen bzw. Scherungen, die aber nicht weiter vorkommen, 
ferner U, V, W die cartesischen Koordinaten des Verschiebungsvektors sowie EZ den Elasti- 
zitätsmodul, # die Querkontraktionszahl und @ den Schubmodul. Hier und im folgenden sind 
dabei diejenigen Terme durch Unterstreichen kenntlich gemacht, die in der klassischen 
Theorie vernachlässigt werden. Der Grund für die Einklammerung von (5c) wird gleich zur 
Sprache kommen. 

Die aus der Balkenbiegungstheorie übernommene einfachste Annahme?) für die Biege- 
spannungen ist ihre lineare Verteilung über die Plattendicke: 


32 32 3% 
0, = M,(@, Na = Mey) zn a N Yo 
die Querschubspannungen sind dann nach einem parabolischen Gesetz verteilt: 


ser, (ey): a _ = i 7 =): af _ eu EH (ZUB) 


Für die senkrechte Normalspannung o, (Randbedingungen: 0, = — p für z= — hund o,— 0 


für z= + h) ergibt sich 
3 RE REN 
u anatrala leeren er 


Von der Verträglichkeit der Ansätze (6), (7), (8) mit (1), (2), (3) überzeugt man sich unmittelbar; 
bei Heranziehung von (4) erhält man natürlich die Gleichgewichtsbedingungen für die Span- 
nungen, die wir hier nicht angeschrieben haben. - 

Die Annahmen (6) über die Biegespannungen machen — wie in der klassischen Theorie — 
die Gleichung (5c) unbrauchbar®), und sie scheidet daher für die weitere Behandlung aus. 

Die Annahmen (6) zusammen mit den Gleichungen (5) nach Fortlassung der unterstrichenen 
Terme sind die formelmäßigen Ausdrücke für die bekannten Näherungsannahmen der Kirch- 
hoffschen Plattentheorie: 

1. Die Mittelfläche der Platte bleibt unausgedehnt (‚‚neutrale Fläche‘). 

2. Erweiterte Bernoullische Annahme: Jede auf der neutralen Ebene vor der Form- 
änderung senkrecht stehende. Gerade bleibt auch nach der Formänderung eine senk- 
rechte Gerade auf der durchgebogenen Fläche. 

3. Der Einfluß der senkrechten Spannungskomponenten auf die Formänderung wird 
vernachlässigt (vgl. S. 162). | 

In entsprechender Weise wie vorhin (vgl. S.162) die Spannungen in z-Richtung zu re- 
sultierenden Kräften und Momenten zusammengefaßt wurden, bilden wir Mittelwerte aus 
den Verschiebungen Z, V, W, indem wir die Summe der Arbeiten der Spannungen 


in den Schnittebenen = const. und y=const. unseres Plattenelementes (vgl. Bild 1) der 


Arbeit der entsprechenden Schnittkräfte und Schnittmomente (vgl. Bild 2) gleich- 
setzen. Dann ergibt sich für eine Schnittebene & = const.: 


hr h 

Sn. We=Q,:w, |, U&k=M,-o,, fü, Y”’d&=H-o, an darB,c). 

—h —h : —n 
®) Vgl. 8. 162. fi 

- 4). Denn diese Annahmen stehen im Widerspruch zu den Kompatibilitätsbedingungen, d.h. streng be- 

trachtet ist das angenommene Spannungssystem in einem isotropen festen Körper gar nicht möglich. Da 


jedoch dieses Spannungssystem die mittlere Spannung über die Dicke richtig darstellt, so wird für die tat- 
_ sächliche Ennnihg ae Annäherung geleistet, falls nur die Platte hinreichend dünn ist. 


11* 


2 R . . = . . Math. Mech. 
164 Schäfer, Über eine Verfeinerung der klassischen Theorie dünner Platten Tue 
BIS ChEN FE NPOTIE EU ir ESS BE SENPeJEBEE 


und analoge Beziehungen hat man für eine Schnittebene y—= const. Durch Einsetzen von (7a) 
in (9a), (6a) in (9b) und (6c) in (9c) folgt 


h h h 
3 '/2\? EB [ 2 in, f z 
ar, | u-i)je a-om Hr (10.227 


— Z 


Die Deutung der Integralbeziehungen (10) ergibt: wist der über die Plattendicke gebildete 
Mittelwert der Verschiebungen W aller Punkte, die auf einer Normalen zur neutralen Fläche 
liegen, und die Winkelgrößen w, und o, sind äquivalent, aber nicht identisch, mit den Neigungs- 
tangenten der Normalen zur neutralen Fläche gegenüber ihrer ursprünglichen Richtung. 


3. Die Grundgleiehungen und die Randbedingungen 

Die noch fehlenden Bestimmungsgleichungen gewinnt?) man nun leicht, indem man die 
Beziehungen (5) zwischen den Verschiebungsänderungen und Spannungen umformt in Be- 
ziehungen zwischen den nach (10) gemittelten Verschiebungen und den Schnittkräften bzw. . 
Schnittmomenten. Zu diesem Zweck werden in (5a), (öb) und (5d) die Ausdrücke (6) und (8) 


eingesetzt, sodann wird beiderseits mit Sn multipliziert und über z integriert, was unter Be- 
rücksichtigung von (10b,c) und Einführung der ‚‚Biegesteifigkeit‘‘ der Platte durch 
3 aa 2 BR 
= 31— en FE a er DE (11) 
[nach entsprechender Zusammenfassung von (5a) und (5b)] auf die drei Gleichungen 
|, „9 , 3nl + 
M= N tu uhr re Re 
,f00y , 20, , 3ull+ a) | 
m—n| dy +4ZZ DE ENmE 3 ee FAssee (12b), 
} _1—u ,„,[00, ei 
H=— se a ee . (12e) 


führt. Zwei weitere Gleichungen folgen aus (5e) und (5f) nach Einführung von (7), beider- 
AR BEE ENRS 3 as 
seitiger Multiplikation von — ji — (z) } wegen (10a) und Integration über z. Beachtet man 


Ah h 
dabei, daß mittels Produktintegration 
h 
Bo 2\? 3 [ z 4 
so ergibt sich 
FREE FA 
0, = — Fries a 7 7 Die. BR ae (12d), 
2, 6 1+u 
wo, = ur + 5 . nd SO Eee We oe, (12e). 


„Aus ren: RER (12) zusammen mit den drei Gleichgewichtsbedin- 
gungen assen sich die acht unbekannten Größen M,, M,,H ß 
ermitteln. . nen ns I N 

. Demgegenüber stellt in der klassischen Theorie (12) (Fortlassen der unterstrichenen 
Glieder!) nicht mehr ein System von fünf unabhängigen Gleichungen dar, da zwischen den _ 


5 : . [N dv, a 
beiden letzten die Beziehung Er Fr besteht. Dieser Sachverhalt bedeutet einen wesent- 


lichen Unterschied zwischen beiden Theorien, dessen Folgen bald klar zutage treten werden. 


Die Elimination von ®, und w, aus (12a,b,c) führt mit Berücksichti 
wegen (11) auf die Gleichungen Re ee 


w 02 w 4 80 2 | 
Meer un NE ne ER nn a A 
£ ar 59 Musi pe 2... (ie), 
a AHEEEERTERE GE HEBEN NERRREERE 
°) E.Reißner [7] leitet diese Gleichungen mittels des Castigliano inzi Der hi 
folgte direkte Weg ist grundsätzlich der gleiche wie ihn A. RE. red = ee 


von verwandten Größen zu einem komplexen Ausdruck — ji i i i 
trägliche Kenntnis der Vf. Herrn W. 2 erna, Bee er ee as ER 
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7 6% 2 are (4 55 h 57 a ee elton), 
0 w 00) 00) 
H=—(1— nee | 
( atslate h SE 


und Einsetzen von (13a) bzw. (13b) und (13c) in (4b) bzw. (4c) mit Berücksichtigung von (4a), 


worauf die Operation Er bzw angewandt wird, führt auf die Gleichungen 


HP 
2 oAw 2 0» 
Sue i 
%— 5 AQn ch Nr ae TR (13d), 
De yon 2 OP ‚s 
Bar Eier en ee öde) 


02 02 
(de satz) 


in denen man für A=0 die Grundgleichungen der klassischen Plattentheorie®) wiedererkennt. 

Aus diesen Gleichungen gewinnt man, indem man die abhängigen Veränderlichen trennt, 
einmal (wie in der klassischen Theorie) eine Differentialgleichung 4. Ordnung für w allein, außer- 
dem aber — und dies ist das entscheidend Neue in der Reißnerschen Theorie —eine weitere 
Differentialgleichung 2. Ordnung für eine Momentfunktion, deren Ableitungen wesentliche 
Zusatzglieder für die Querkräfte Q, und Q, liefern. 

Um die erste grundlegende Gleichung für w zu erhalten, setzen wir (13d) und (13e) in (4a) 
ein und berücksichtigen dabei den durch Anwendung des Operators A auf (4a) gegebenen Zu- 
sammenhang. Die sich ergebende Bipotentialgleichung 

i 


N A4lu=p— Eh An: ae RAR TEE (D 
sta) 


unterscheidet sich formal von der entsprechenden in der klassischen Theorie nur durch einen 
unwesentlichen Summanden des inhomogenen Gliedes. Freilich ist die Bedeutung von w nicht 
in beiden Fällen die gleiche: In der klassischen Theorie bedeutet w die Durchbiegung der neu- 
tralen Fläche, in der verfeinerten Theorie dagegen den Mittelwert der Vertikalverschiebungen 
längs einer Normalen zur Mittelfläche (vgl. S.164). 


Setzt man | 
EN RU RR Te (14), 

wobei w* eine Partikularlösung von (I) sei: 
| Nam, Ap- a ae es): 


so ist ® eine willkürliche biharmonische Funktion: 


[ 4148=0 | a N ee (1). 


Ferner suchen wir die Gleichgewichtsbedingung a durch den Ansatz 


„?4P 
Q; — or nn nr u? n AS IH TER SE NE (15 a), 
3 
Q,= 0% ZN: 2 - Fr BR Er gen Meer (15b) 


unter Einführung der Momentfunktion y zu befriedigen, wobei Q und Q, Partikularlösungen 
der Gleichungen i 


Bey .e dm" 2 2% Ps FERSAIIERE 

bzw. 3 
92 2 En 9Aw* 20 89, ar 13e* 
ren 


6) Vgl.\etwa [9]; s. 212. ‘Formeln (5) und (7). 
sa Er‘ | 


3 : = Z.angew. Math. Mech. 
166 Schäfer, Über eine Verfeinerung der klassischen Theorie dünner Platten 34.92 Nr. 6 Juni 1952 
ee tee se 
RN ee en Er 


sein sollen. Differenziert man (13d*) partiell nach x und (13e*) nach y und addiert beide Glei- 
chungen, so folgt mit (I*) in symbolischer Schreibweise 


1-30 End 


0% 0y 

was durch & 
N TR (16) 
0% 0y 


ü i ') wird i befriedigt. 
erfüllt wird. Wegen (16) und (I’) wird in der Tat (4a) durch den Ansatz (15) 
Einsetzen von (15a) in (13d) und (15b) in (13e) liefert dann wegen (14), (T‘) und (134*) 
bzw. (13e*) i 
0 2 15 \- 214) =0 
sr) ern) 0: 
woraus folgt, daß der Klammerausdruck eine Konstante sein muß, die gleich Null gesetzt werden 
kann. Daher ergibt sich als zweite grundlegende Gleichung: 


5 
A—zayal| ArbN. on. nenn (II). 


Hinsichtlich der eindeutigen Lösbarkeit der Randwertaufgabe dieser elliptischen Differential- 


gleichung ist entscheidend, daß der Faktor ——, eine positive Zahl ist”). 


2.h? 
e dx Für die beiden Differentialgleichungen (I‘) und (II) können drei 
3 dy VEr@ unabhängige Randbedingungen gestellt werden. Indem wir auf die - 


Darstellung (12) zurückgreifen, können also am Rand der Platte (vgl. 

Bild 3; an Stelle der cartesischen x-y-Koordinaten sind hier natürliche 

Koordinaten in Richtung der Bogenlänge (s) und der äußeren Nor- 

4) malen (n) der Randkurve gewählt) auf Grund der verfeinerten Theorie 
.  die.drei Bedingungen: 


| Uns. u a eh 
oder 

Bild 3 Q,=9,,;, M,.=M,, H=H,‘... - (17b) 
vorgeschrieben werden; so z.B. für die eingespannte Platte: v=&,=@, —=0 oder für den 
freien Plattenrand: Q,—= M„=H,—=0, während bei der frei gelagerten Platte »— 0 an die 

Stelle von Q, = 0 tritt. i | 
In der klassischen Theorie artet (II) in die nichtssagende Beziehung y = 0 aus, das Inte- 
grationsproblem ist um zwei Ordnungen erniedrigt, und daher kann man eben dort — unter 
Umständen im Widerspruch zu den Grundforderungen der Elastostatik — nur zwei unab- 
hängige Randbedingungen befriedigen. Wie dieser heikle Punkt durch das Theorem von Thom- 


son und Tait in zwar sinnreicher, aber doch etwas künstlicher Weise überbrückt wird, kann 
hier unter Verweis auf das einschlägige Schrifttum®) unerörtert bleiben. 
! 


4. Anwendung auf die frei gelagerte Rechteckplatte 


Als Anwendung der Reißnerschen Theorie wollen wir eine auf Biegung beanspruchte _ 
Rechteckplatte (Kantenlängen 2@ und 25) untersuchen, die an den Rändern frei gestützt wird 
und nach dem doppelten Kosinusgesetz | REEL 


r 


pP= m C0S 3 cos =, SE aa ee ee wahre (18) 
belastet ist (vgl. Bild 4, wo im Grundriß die Linien gleicher Belastung eingetragen sind). 
Die Randbedingungen unseres Problems (vgl. Bild 2) sind®): j ’ 
v=M, =-H=0r: fr, 2er See (19a), 
v=M,=H=0 für. ya re .. (19b). 


?) Vgl. etwa [10], Kp. IV, $5,1. x 
8) ur u.a. BL S. 212f. oder auch [11], 8. 78 £f. 
a 


; ir haben dabei naturgemäß längs des Randes das Verschwinden der gemittelten Vertikalver- 
schiebungen zu fordern. £ 
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Wir beschränken uns an dieser Stelle hauptsächlich auf die Durchführung der Rechnung 
und verzichten auf eine eingehendere Diskussion unseres Ergebnisses, insbesondere auf einen 
quantitativen Vergleich dieser verbesserten Lösung 
mit der sich aus der Kirchhoffschen Theorie er- 
gebenden!?), 1 EBR 

Als mit den Randbedingungen verträgliche Parti- A DR 
kularlösung von (I) setzen wir im Hinblick auf (18) an: eek 


EV 
> cos Be (20), 


w* = Wy COS 


wonach sich die Konstante w, aus (I*) mit (18) be- 
stimmt: 


x 16 a® bt 


=» 1 | ( u) (a? ) pe 
| 5(l—u) 4a?b2 


Bild 4 


In entsprechender Weise werden Partikularlösungen für die Querkräfte angesetzt. 
EEE FE LER 
De sin 3, 08 >73 


eg, cos sin U | 
0, = 9 cos a a EN (22b). 
Einsetzen von (22a) in (13d*) und (22b) in (13e*) ergibt 


zart b) „u (a? -+ b2) 2 RER 
alıtz u an a Bee en 


bzw. 


2% la — b2 3(g2 4 b2 9 BERAR 
En ED een ee Ne 


‚woraus gemäß (16) wegen (21) 


folgt. 

Unsere Partikularlösungen (20), (22) entsprechen der endgültigen Lösung in der klassischen 
Theorie. Hinsichtlich der Durchbjiegung w und des Biegemomentes M erfüllen sie bereits die 
‘ Randbedingung, nicht aber für das Drillmoment Z. In diesem Punkt greift nun die ver- 
feinerte Theorie entscheidend ein. 

Zur Befriedigung von (I’) mit (19) bieten sich wegen der doppelten Symmetrie des Problems 
als Partikularlösungen die biharmonischen Funktionen 

vru|, go "UL Br Sin ?73 

TER ER: ı EE K ER RAR 24) 

vi VE VNE. VC% ( 
+ cos = a,&oj Dr + b, 9b Sin = 


AT er 


an. . ® . 
Im Hinblick auf (15) mit (22) müssen die Partikularlösungen, die wir zur Befriedigung von 
(II) mit (19) ansetzen, doppelt ungerade Funktionen sein: 


2 U; EYRY 
| w=(, a er Sin A,y+ 0, sin Er I AS ee er Vie (25), 
wobei = Er 
vn 5 fe, 5 = 95%) 
= A,= 1a t3% und Ne P Am ram 4, > Be) ( 
gilt, : 


10) Vgl. hierzu etwa [11], 8. 73ff. 


& ENT ERETR 
: ar 
SE | 
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Als Lösungen von (I‘) bzw. (II) setzen wir zunächst rein formal die unendlichen Reihen 


d— Di, 2 SR ee a (26) 
y=1,3,5, 
und 
ya DE we ee ee (27) 
y=1,3,5,... 


an. Die Freiwerte A,, a,, B,, b,, und O,, c, dieser Entwicklungen werden so bestimmt, daß die 
Randbedingungen (19) erfüllt sind. 
Zu diesem Zweck denke man sich (24) in (26) und (25) in (27) eingesetzt. Die beiden Rand- 
bedingungen w= 0 und M=0 zufolge (13a) bzw. (13b) mit (14), (20) und (15a), (22a) bzw. 
92a * 271% * * 
(15b), (22b) sowie (18) Fr ; = ‚p und = bzw. = verschwinden dabei von vornherein 
-am Rand) liefern dann Cosinus-Reihen, die identisch verschwinden sollen. Folglich muß jeder 


einzelne Reihenkoeffizient Null sein, und das ergibt für jede ungerade Zahl » die zwei 
Gleichungspaare: 


v=l für = +.e: weh Ü on: Henne (28a), 


GEIT: 
v=0 für y-+b: a a (28b); 
4 a ren ; 


MN: Tür’ 2-0; 


vanı? Va va | i 
0 N la, (1 -u)Cof te (2 + Sin n 
a 
8 y2rıı VITA 4 VIE ’ 
Ta De s »)}+o, - a] ),a* h? 
MN für y=-+b: i 
2 r2 
Te Se vb ( vb vab_.. vb 
Tan [Ar UWE + B [260] Zr + u) Ein x 
8 vn? vb ) 4 vn 
51 &l 4b, 557 A,Cof A,b- he 


Dagegen liefert die dritte Randbedingung H = 0 zutage (13c) mit (14), (20 und 15a 
(22a) bzw. (15b), (22b) vorerst nur ein Paar von Summengleichungen, ai ar ER 
Konvergenz der Reihenentwicklungen vorausgesetzt — die Sinus-Glieder auf die linke Seite 
abgesondert seien: u 


H—=09-fürtz= a: 


> Ja-WnE +0, Ein ?7E 15, (Sin?22 1 PRRg, Pra. 
an a BIER ah 


VE UNE UN TR N 
EreenrUnerz vl 5 + ge) Ein Aa > 
2, M\.zell.sin ?*Y | | RER 
seht) || RZRR (30a), 
ns 2° (u- ne |4,.Cin?”V RYLE 
Et Na [Ar Ein + Bu (amt Yo | 


I yry‘ a 
MENT Arne »)| +0, (6 +) ein Ay: 1) 


vr 
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= 2 772 
ee | ‚vnb vnb..vnb 
>, . 4.8 Wr + 4, Sin Z, + B, -\&in > ' = Cof 5 


2178, Done 
8 ven? _. vb 2 27,2\ 
SnTeRrr er nl (a4) NE 
DER VIE 
z n(55+ TE 207 (0b): 


© vrl 272 
nn 112 ka- ver | en ( VIE VIER VIE 
| ) MN Tr |® >» Su +, 


8 var? _, VIL 2 y2gt2\ _, 
era a ee: 3% a) Sin Ava hr} 


bei der Summation ist zu beachten, daß die eingerahmten Glieder 


272 
ER IaTz war SE 
( ENT pe 5 2b" 2a ” 


nur für v=1 auftreten, während sie für v» >1 fortzulassen sind. 


u Um nun auch die aus (30) fließenden Beziehungen zwischen den Freiwerten ablesen zu 
ash müssen zuvor noch die rechten Seiten der Gleichungen (30a) und (30b) in Fouriersche 
eihen entwickelt werden. Man hat dafür folgende Entwicklungen heranzuziehen, wobei wir 


vorübergeh i RT T EN b 
ergehend die Bezeichnungen 1a > Pr Tr ed? Wem hY— Ne 


A,@a=Gyı, A,b= &,s einführen: 


NE 
A “ n—ı sinn 
Sin v&,—= —va,Co| En >: (— 1) Be 
id ad, n? + v? 0% 
; in ER 
4 < "Z) m—r?of TC, 3 
ve,&ojv&,—= — vo, —1) 2 Eee ee ee Z 
| 2 22 Er (n?-+v20,)? Pan 3 na (31) 
Rn: 
, 8 r 6) Pal sın ER 
Sinn, = — %,, 00 eh re 
Me 103 en a 40, 
28,Dyeo» nt ” 
12 
(#1, 2% ee ER ER! 


Durch Kombination von (28a) und (29a) bzw. (28b) und (29b) ergibt sich 


2N ca IE var? vn ya var i 
lt te ne 
Ein Zz R 
= 0:54, 6oj A,a:h? 
ler N ara P Blue vb 1 van? 12 
Es emziaen j=Bs Ar! re (32b). 


= 0,:5 4,60) A,b- he 


Ehe wir an die Bestimmung der c, und (, aus (30) gehen, können wir uns auf Grund der 
ersten Voraussetzung (h< a,b) eine Vereinfachung der Ausdrücke (32) erlauben. Indem wir 
zunächst einmal große Werte von » außer Betracht lassen, dürfen wir für alle v, die kleiner als 
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ein gewisses », sind, näherungsweise v?h? gegen a? und 52 unterdrücken. Daher dürfen wir in 


(25*) für »<v, setzen: 
Ash area N RS ec (33), 


ferner an Stelle von (32): 


VITA VITA VI ya 
ren. ae (4a) 
a aCofAa 12° Cofic h : 
N N 
ER a Van 2. ..(34b) 
Ara 22 aCojAa EB TEL NIEREN 


Hat man die Entwicklungen (31) in (30) eingeführt, so ergibt der Koeffizientenvergleich 
entsprechender Sinus-Glieder ein unendliches Gleichungssystem für die o,, O,, das zwar grund- 
sätzlich nicht rekursiv auflösbar ist, sich aber deswegen leicht iterativ behandeln läßt, weil in 
jeder Gleichung zwei betonte Glieder auftreten (d. h. Glieder, die sich durch ihre Größenordnung 
stark von den übrigen abheben), welche näherungsweise einander gleichgesetzt werden können. 


Wir.führen den Koeffizientenvergleich für sin 3 und sin >, 
schenrechnung erhält man = 


5Al—-u)a®Nw hr 
2 4ab 


| - 2b er 
— fa —u) = -aGofaa-zg Te 4 — EofAa ; 395, E= acofaa]— Vesinaanal.c 


durch. Nach kurzer Zwi- 


(35a) 


< —= CofAb | 85A2CofAb 
Be 172 fa (d—u) ba aaa 7232 Fe "Cr 
bzw. 
5( — u)n?Nw ja 
2 4ab 


2 - 5 

— fa ga |560f Ab Tg? + — Cof 1-29 3° +5Cof + zEinadan).c, (35b). 
— \ = CojfAa | 8aRrkoj Aa 

FE Er Et] 


Da v3 hY=% sehr groß ist, folgt sofort 


(1— u)? Nw, ; 
UT abe ia Ne Sa 
(1 — u) n2Nw, 
ATI a re Ser 


und hieraus weiter nach (34a) bzw. (34b) und wegen TgAa =TgAb 1: 
SH a ie 
tw, 


BEE | 
ee PET Ze N ee > BE uch N! 
(1-u)m2Xg = 
2 3 DR 
4-22 Bl -y3: Dun DE (37b) 
4aßoj 2bCoj 5 ä 


Die weiteren Gleichungspaare brauchen wir gar nicht explizi i ie si 5 
\ plizit anzuschreiben. Sie sind ent- 
sprechend (35) aufgebaut, wobei aber nun wesentlich ist, daß in der »-ten Gleichung (Koeffizienten- 


vergleich von sin v&,(@=1,2;9=3,5,,...)statt co, in (35a) c, und statt CO, in (35b) C, steht und 


die linken Seiten Null sind, weil für » > 1 in (30 die I-Gli 
SE (30) die Z1-Glieder fortfallen. Man erkennt dann 


daß alle übrigen Freiwerte im Vergleich zu den Freiwerten c, und Cjeine 
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Größenordnung 7 kleiner sind. Wir haben deshalb auf die Wiedergabe weiterer Bestimmungs- 
gleichungen verzichtet, die aus dem Koeffizientenvergleich der folgenden Sinusglieder resul- 
tieren; denn es sei dahingestellt, ob in einer linearen Theorie (um die es sich im vorliegenden 
Fall ja handelt) die Entwicklung. der Freiwerte, die — wie eine genauere Untersuchung lehrt — 
Potenzreihen in Ah sind, über die 1. Potenz hinaus überhaupt sinnvoll ist. 


5. Ergebnis und Diskussion 


Tragen wir die gefundenen Ausdrücke (36), (37) für unsere Freiwerte in (24) bzw. (25) ein, 
so erhalten wir mit (14), (20), (21) bzw. mit (15), (21), (22), (23) als Lösung unseres Problems 


| no ein” 


LEE EL 20, y__2a| „me 
EEE iR @ re, nl 2a 
Cor > Co” 
. 2a 2a 
Ta TT% 7% = i 
[2055 &o1 >73 Ein 577 2 (38 
fr - — cos X |, 
b & Ta ads Be 2b 
055 2b 
Tat N 
und 
Gone nY 2(1—wab{1/5 . nabojiy b m my Sin Ae| 
Ge aa ala 9) 2" 2a Sina 2° Einza|’ er 
a 
3 2 02b 
19T z(@® + 2) P° 
bzw 
Q a2. ay 2ll—mabli/5 . nylolia a m a« Sin Ay] 
en ne HB) IV 2° 2 india h 2 "2a Sin äbj’ oh 


> 2ab? 
HATT? = 


Deuten wir kurz dieses Ergebnis, so zeigt sich einmal, daß die verfeinerte Theorie auf die 
Durchbiegung der Platte keinen bedeutenden Einfluß hat, was nicht überraschend EN 
Dagegen hat die Querkraft eine wesentliche Modifikation erfahren, die wegen des Faktors 1/ 
besonders in der Umgebung der Ecken entscheidendes Gewicht bekommt und dort eine 
erhebliche Anpreßkraft ergibt. Bekanntlich versagt ja die klassische Lösung hier. Ein starkes 
Anwachsen der Querkräfte in den Ecken ist übrigens bereits von Experimenten her bekannt. 

Abschließend möchten wir sagen, daß die Reißnersche Theorie eine — wie uns scheint — 
recht fruchtbare Erweiterung und Verfeinerung der Kirchhoffschen Theorie darstellt, deren 
weitere Bewährung man mit Interesse verfolgen wird. Auf jeden Fall bietet diese neue Sa 
den großen Vorteil, daß sie in natürlicher Weise die Erfüllung aller erforderlichen dm: 3 
bedingungen ermöglicht und hierfür keiner zusätzlichen Hilfsbetrachtung bedarf, wie dies die 
klassische Theorie nötig hat. 
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Grenzschichten an zylindrischen Körpern mit nichtstationärer 
Querbewegung 


Von Walter :Wuest in Göttingen 


Ein besonders einfacher Fall von nichtstationären dreidimensionalen Strömungsvorgängen in laminaren 
Grenzschichten ergibt sich bei der Anströmung zylindrischer Körper mit überlagerter zeitabhängiger Quer- 
bewegung. Ein allgemeines Lösungsverfahren wird angegeben und an mehreren numerisch berechneten 
Beispielen erläutert. 3 

A very simple case of threedimensional nonstationary problems of laminar boundary layer currents 
is given by the flow against cylindric bodies with axial movement depending from the time. A general method 
of solution is deduced.and discussed for several numerical examples. 


Un cas tre&s simple de probl&mes d’&coulements troisdimensionels nonstationnaires dans des couches 
limites laminaires se presente par l’&coulement autour de corps cylindriques soumis & un mouvement axial 
dependant du temps. Une methode generale de solution est donnde et discute pour plusieurs exemples 
calcules numeriquement. 

OcodeHHo UPocToä CIyYali HeYCTAHOBHBIIUXCH TPeXMePHEIX Pe3KHMOB IIOTOKA B JIAMUHAPHBIX 
TOTPaHHYHEIX CHOAX TOAYAYAETCA IIPH HATeKAHHHY HA IIHIIHHNPHYeCKHE Tea C HAC/IOCHHBIM HONe- 
PeUHEIM ABH>KeHHEM, BABHCHMEIM OT BpeMenHn. Jlaerca od METOA peIleHus, IOACHAeMBHÄ HA 
HECKOJIBKHX YHCJIEHHO PeINeHHBIX IIPHMepaXx. : 


1. Einführung 

Bei der rechnerischen Untersuchung von laminaren Grenzschichten hat man sich bisher 
h fast ausschließlich auf ebene Strömungsprobleme beschränkt, bei denen innerhalb der Grenz- 
schicht nur zwei Geschwindigkeitskomponenten auftreten. Soweit man eine zeitliche Veränder- 
lichkeit in Betracht gezogen hat, wurden im wesentlichen Anlaufströmungen behandelt, bei 
denen ein Körper, beispielsweise eine ebene Platte oder ein elliptischer Zylinder, aus der Ruhe 
heraus in Bewegung gesetzt wird [1], [2]. Diese Untersuchungen sind auch auf den, Fall aus- 
gedehnt worden, daß ein Körper aus einem bestimmten Bewegungszustand heraus beschleunigt 
oder verzögert wird [3]. 


Bei Grenzschichten in dreidimensionalen Strömungen treten innerhalb der Grenzschicht 
7 drei Geschwindigkeitskomponenten in Erscheinung. Wegen der rechnerischen Schwierigkeiten 
hat man sich bisher auf einige wenige einfache Fälle beschränkt. Eine besondere Vereinfachung 
r tritt dann ein, wenn die drei Geschwindigkeitskomponenten nur von zwei Koordinaten abhängen. 
- Derartige Aufgaben treten bei der schrägen Anströmung von zylindrischen Körpern auf [4]. Läßt 
3 man etwa die Geschwindigkeitskomponente v in die Richtung der Zylinderachse fallen, so ist 
Nr die Geschwindigkeit u, die in Richtung der Mantelfläche senkrecht zu v verläuft, unabhängig 
"iR .von vund kann nach bekannten Verfahren der zweidimensionalen Grenzschichttheorie ermittelt 
E67 werden. Hinterher kann man v aus einer linearen partiellen Differentialgleichung ermitteln, die 
(a in gewissen Fällen auf eine gewöhnliche Differentialgleichung zurückgeführt werden kann. Für 
BD ‘ die Staupunktsströmung hat L. Prandtl [4] die Lösung berechnet. 
ne _ In Erweiterung der Prandtlschen Untersuchungen sollen im folgenden Strömungsvor- 
Br. gänge an zylindrischen Körpern untersucht werden, bei denen diese in Richtung der Zylinder- 
Be achse zeitlich veränderliche Bewegungen ausführen. Neben der Staupunktsströmung werden 
; insbesondere Strömungen .an keilförmigen Körpern und am Kreiszylinder behandelt. 


| - 2. Grenzschichtgleichungen und allgemeine Lösungen 
Be . Wir betrachten die Strömung um zylindrische Körper, die sich beiderseits ins Unendliche 
BR erstrecken. Die abgewickelte Mantelfläche bilde die &*, y*-Ebene, wobei die y*-Richtung mit 
der Zylinderachse gleichlaufend ist und x*=0 mit der Vorderkante 
(Staulinie) zusammenfällt. Die z*-Richtung verlaufe senkrecht zur Wand. 
Der Körper werde entsprechend Bild 1 senkrecht angeströmt, voll- 
führe aber gleichzeitig in Richtung der Zylinderachse Bewegungen 
mit zeitlich veränderlicher Geschwindigkeit. An der Wand (d.h. für 
2*=0) ist daher = w*—=0,v*=YV*(t), am Außenrand der 
Grenzschicht dagegen sei u*—= U*(a), v*—=0. Die Geschwindigkeits- 
komponenten u*, v*, w* sind offenbar von 9* nicht abhängig, so 
un: Btegnunng Ben daß die Ableitungen nach %* in den Grenzschichtgleichungen ver- 
schwinden. Durch Einführung der Anströmgeschwindigkeit U undeiner 


‚passend gewählten festen Länge Z kann die Reynoldssche Zahl R=U, ei gebildet 
v 


werden und die Grenzschichtgleichungen können durch folgende Umformungen dimensions- 


! » 
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los und unabhängig von der Re-Zahl gemacht werden: 


ur v* EN 
ee) Re 
DE U 
ee Ey 
U D> ti=t ‚ 
* zr\, — 
==. — — I_ / 
& I 2 (Z)/Re. 


Damit gehen die Grenzschichtgleichungen über in: 
ou ou ou dU Hu 


I) 0v ov Ov 
& a 9% Ar 92 022 2) 
ou  ow 
ee u en Be ao 
Man erkennt zunächst, daß die erste Differentialgleichung v nicht enthält, so. daß kein Grund 
für eine zeitliche Abhängigkeit von u besteht. Daher kann = —=(0 gesetzt werden. Wegen 


der Kontinuitätsgleichung (3)ist auch wnicht abhängig von der Zeit. Die Gleichungen (1) und (3) 
sind mit den bekannten Grenzschichtgleichungen für den ebenen stationären Fall identisch und 
können demnach nach bekannten Verfahren gelöst werden. Da somit u und w vorab ermittelt 
werden können, ist Gl.(2) eine lineare partielle Differentialgleichung für. v(z, 2, t). Die Ab- 
hängigkeit von der Zeit tläßt sich abtrennen, wenn man für die Querbewegung folgenden Ansatz 
machen kann: 
Mr V,er% 

sowie für die Grenzschichtströmung 

Der 9m teuer N zand): 


Wegen der Linearität der Differentialgleichung (2) kann man solche Lösungen ohne weiteres 
überlagern, so daß auch 


v5 Vyentr, (2,2). Es ne (4a) 
% F 


eine Lösung darstellt. 

Ein wichtiger Sonderfall ist die Grenzschichtströmung an einem periodisch schwingenden 
zylindrischen Körper. Bei Beschränkung auf eine einfache Sinus-Schwingung setzt man 9,— +10 
und für die Querbewegung 


7 (eiot 4 0-0) = Vycosot Rn ee er N) 


und für die Grenzschichtströmung 

=, (5) 0itn(n)sinoi). .... vn. 98) 
mit 

4 &0)=1; ,(,0)=0; (a, @)=7(m, @)=0, 
Bei ebenen Grenzschichtströmungen erhält man unter gewissen Voraussetzungen ‚affine Grenz- 
schichtprofile‘‘, bei denen die Abhängigkeit von-der «-Koordinate nur als Maßstabsfaktor in 
Erscheinung tritt, während die Grenzschichtgleichung auf eine gewöhnliche Differentialgleichung 
zurü:kgeführt werden kann. Es erhebt sich die Frage, ob es auch im vorliegenden Fall affine 
Lö ungen für die v-Komponente gibt, bei denen die Abhängigkeit von x nur in einem Maß- 
stabsfaktor m(x) zum Ausdruck kommt, so daß also —r(mz). Wenn wir die Lösung der 
Gleichung (2) in der Form > 
; v— V,e* r(mz) 

ansetzen, erhalten wir: 
arme) um! umi=mt! ..... 2er (6). 
_ Da «eine Konstante ist, müssen auch um’ und wm sowie m? Konstante sein. Einen konstanten 
Maßstabsfaktor m erhält man aber nur bei der Staupunktsströmung. In allen anderen Fällen 
ergeben sich aber auch bei affinem «-Profil keine affinen v-Profile. 


= x STREBEN, 
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Folgende Fälle werden näher untersucht: 

1. Staupunktsströmung (u-Profil und v-Profil von & unabhängig). 

2. Keilströmung (u-Profil affin, v-Profil nicht affin). Sonderfall: Längsangeströmte ebene 
Platte, die in ihrer Ebene aber senkrecht zur Strömungsrichtung instationär bewegt wird. 

3. Allgemeiner Fall (weder u-Profil, noch v-Profil affin). Sonderfall: Angeströmter in 
Achsenrichtung instationär bewegter Kreiszylinder. 


3. Staupunktsströmung 

Die Staupunktsströmung bildet sich an einer quer zur Anströmrichtung gestellten Platte 
aus, die in der einen Richtung unendlich lang ist. Allgemein tritt die Staupunktsströmung in 
a der Umgebung der Staulinie von 
G=Vez =VeRe Z"/L zylindrischen vorn abgerundeten 
Körpern auf. Für die Geschwin- 
digkeit U außerhalb der Grenz- 
schicht setzt man hier U=ca. 
Die «-Komponente der Grenz- 
schichtströmung ist dann in be- 

kannter Weise gegeben durch: 


u=Uo' (Vez) ; 
wobei o’ eine bei Hartree [5] ta- 
bellierte Funktion ist. Aus der 
Kontinuitätsgleichung (3) folgt 
dann für w 


ER 


so daß die Differentialgleichung (6) 
übergeht in: 


tor ı=0 a 


wobei die Striche Ableitungen nach 
£=Yez bedeuten. Gl. (7) wird am 
besten numerisch integriert wobei, 
zunächst o aus der von Hartree 
tabellierten - Funktion 0’ durch 
Integration gewonnen wird. Die 
schrittweise numerische Berech- 
nung der Funktion z(£) ist für die 


Werte — —=1 und = = —-I durch- 


e 
ee ulu,v/v 


0 02 04 06 08 10 12 14 geführt und in Bild 2 aufgetragen. 


ES 
Bild 2. Die Geschwindigkeitskomponenten der Grenzschichtströmung an einer Für — =0 (konstante Querge- 
beschleunigten oder verzögerten Platte (Staupunktsströmung } c 


schwindigkeit ist die Kurve der 
Arbeit von L. Prandtl [4] entnommen worden. Das Grenzschichtprofil der u-Komponente: 
ist nach Hartree [5] eingetragen. 


Wegen der linearen Überlagerungsfähigkeit der Lösungen für v kann man auch sofort 
den Fall V=V,+ Vje’'+-Vze”°t behandeln. Insbesondere können die berechneten Lösungen 
auf folgende Fälle übertragen werden: V=V,(1— et); V=Y, Sinct und V=V,6of ct. 

Die Rechnung ist ferner für die Staupunktsströmung an einer schwingenden Platte durch- 
geführt worden. Für die Bewegung der Platte ist dabei der Ansatz (5), für die v-Komponente 
der Grenzschichtströmung der Ansatz (5a) zugrunde gelegt. Dabei sind aber die Funktionen 2; 


und 7, nur von {=Jez abhängig und durch folgendes gekoppelte System von gewöhnlichen 


Differentialgleichungen bestimmt: 
da-+ u A 0 | 


toten] 
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Die Sausbedingungen lauten: 
url 1 r,(0) = 0, 
le ER eo) Fe A ER 


Für den Fall 2 Flez 1 sind die Funktionen 7, und r, schrittweise numerisch berechnet und in 


Bild 3 aufgetragen worden. In Bild 4 ist der Gesamtverlauf der Komponente » für verschiedene 


G=YVez =Vere Z/L 


Zeiten aufgetragen. 


5 =Ve'z = VeRez*]L 


TT 
"1,2 Bi 0 1 Vo 
Bild 3. Hilfsfunktionen für die Berechnung der Grenzschicht Bild 4. Die v-Komponente der Grenzschichtströmung an einer 
an einer schwingenden Platte (Staupunktsströmung) schwingenden Platte (Staupunktsströmung) zu verschiedenen 


Zeitpunkten der Schwingungsperiode 


4. Keilströmung 


Bei affinen Lösungen der ebenen Grenzschichtgleichungen hat die u-Komponente bis auf 
Maßstabsfaktoren an jeder Stelle x gleichen Verlauf in Abhängigkeit vom Wandabstand. Man 
kann also setzen: 


u=U(e)o'(m2)=Ul(a)doldE wobei Z=m(a)2.:...... RESAUU: 
Wegen der Kontinuitätsgleichung (3) ist dann die w-Komponente gegeben durch: 
u lo 0% los—to). a RR ER (10). 
dz m m m 
Mit den Ansätzen (9) und (10) geht die Differentialgleichung (1) über in: 
| z Ber Boll al en hear (11) 
wobei e 
lad Se ) Uhr dU 
— mda 2? m? da 


"Konstante sind, so daß also U und. Mm nicht frei wählbar sind, sondern mit der Konstanten c 


er ier sind durch: 


Mn = 2 
Teck) > Aueh), U RR NG (12). 


Die zus seukalelichıng @) für die v-Komponente. geht mit dem Ansatz 


Pu Le (%, ©) 
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x EP mn Tr Te ET 
über in: ER 
re 00. 
a TE g=const ( ) 
Bei Betrachtung der Koeffizienten dieser Differentialgleichung und der Beziehungen (12) ist es 
ad, —(a 


naheliegend, von & zur neuen Veränderlichen & = (# — 2) “® überzugehen und für z(£, &) 
folgenden Ansatz zu machen: 


MOD EERNOE HOT EA ee et) 


Dieser Ansatz führt zu folgendem System von Differentialgleichungen: - 


24, —4 


wtaou—na mi —2nly —W)ern=0.....-- (15), 


wobei a 


Bei der Staupunktsströmung ist a,—@,, also &= const, so daß der Ansatz (14) seinen Sinn 

verliert. Ein ee Sonderfall Ba sich für 2 @,—= a,. Dieser Fall ist nach Mangler [6] 

dann verwirklicht, wenn die Geschwindigkeit zr 

der Potentialströmung am Außenrand der 4 ax Vre & 

Grenzschicht nach einer Exponentialfunktion 

in.x verläuft. & muß also in diesem Fall eben- 

falls als Exponentialfunktion in x angesetzt 

werden und der Ansatz (14) führt zueinem ähn- 

lichen System von Differentialgleichungen wie 
R; (15). Doch soll auf diesen Sonderfall nicht 

i näher eingegangen werden. 

‘Die affinen Grenzschichtprofile, die der 
Differentialgleichung (11) genügen, werden 
durch die symmetrische Anströmung eines Kei- 
les vom Schnittwinkel 8x verwirklicht. Dabei 


x * 
Anströmung 


GR, 


Bild 6. Hilfsfunktion für die Berechnung der Grenzschichtströmung 
an einer Platte (Blasius-Grenzschicht), die seitlich beschleunigt 
bewegt wird 
ist a=1und a,— ß zu setzen (vgl. z.B. Mangler [6])!). ßnimmt dabei Werte zwischen 0 (Platten- 
strömung) und 1 (Staupunktsströmung) an. Der Wertebereich 1< ß<2 entspricht der Strömung 
in einen hohlen Keil hinein. Nach Mangler [6] können andere Wertebereiche von ß als Strö- 
mungen in konvergenten und divergenten Kanälen verwirklicht werden. Die Funktion o’ ist 
bei Hartree [5] für verschiedene Werte von ß tabelliert. Es bereitet daher keine Schwierigkeit, 

das Gleichungssystem (15) numerisch zu integrieren. 

Die Rechnung ist für den Keilwinkel B=0 durchgeführt worden, d.h. für eine längs ange- 
strömte ebene Platte (Blasiusgrenzschicht), die zusätzlich quer zur Strömungsrichtung in ihrer 
Ebene nach einer exponentiellen Zeitfunktion bewegt wird. Die Konstante c in (15) wird dabei 
c=1, wenn man die tatsächliche (nicht dimensionslose) Potentialgeschwindigkeit U*=U setzt. 


Als Rechenbeispiel wurden für den Fall Y= V,e'die Funktionen 7, und r, berechnet und in 
Bild 6 aufgetragen, (or 


5. Strömung an einem Kreiszylinder x er x 
Bei den bisher behandelten Problemen — Staupunktsströmung und Keilströmung — war 
die u-Komponente der Grenzschichtströmung stets affin, also abgesehen vom Geschwindig- 


\) Eine der beiden Konstanten a, oder Ga, kann beliebig angenommen, also auch gleich 1 gesetzt werden, 


Bild}5. Symmetrische Anströmung 
eines Keiles 


h te 
1, - Pe 
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keitsmaßstab nur von £ abhängig. Bei der Strömung um einen Kreiszylinder ist aber bereits 
die u-Komponente nicht mehr affin. Wir kehren daher zu den Ausführungen des 2. Abschnittes 
zurück. Als charakteristische Geschwindigkeit U, wählen wir die ungestörte Anströmgeschwin- 
digkeit des Zylinders und als charakteristische Länge den Durchmesser des Zylinders. Die 
dimensionslose Potentialgeschwindigkeit U am Außenrand der Grenzschicht kann daher in fol- 


gende Reihe entwickelt werden: 
De a (16) 
da die geraden Glieder aus Symmetriegründen wegfallen. Hierbei bezeichnet demnach x die 


vom Staupunkt an gerechnete Abwicklung eines Zylinders vom Durchmesser 1. Die Potential- 
theorie würde ergeben; 


Ur=2sn2T, 
Experimentell hat K. Hiemenz [7] jedoch gefunden: 
NE I R  BES Re (16a), 


wenn man die Ergebnisse von Hiemenz mit einem Zylinderdurchmesser von L= 9,76 cm und 
einer ungestörten Geschwindigkeit U, —=19cm/s dimensionslos macht. Die Formel (16a) 
gilt nach den Versuchen bis x = 0,72. 


Die Geschwindigkeitskomponenten « und w können in folgender Weise in eine Reihe 
entwickelt werden: 


u=aHl)e rs) RH 2... ....(17) 
= —— (110,430, +5 40024: ) ARTE: (18) 
Ya 


mit = Yc2= Je, Re z*/L.-Die Funktionen o, und 0, (bzw. o/ und 0/) sind universell berechnet 
worden. Dagegen hängt der Verlauf von o, vom Verhältnis c,/c, ab?). 
Der nach einer exponentiellen Zeitfunktion bewegte Zylinder. 

Für die v-Komponente machen wir aus Symmetriegründen folgenden Ansatz: 


v= Net SO HLÖAH]) - 0. (19). 
Dies führt zu folgendem Gleichungssystem für die Funktionen z, und 7,, wenn man (19) indie 
Grenzschichtgleichung (2) einsetzt. 


” ’ & 
n+ nu — nz De RETTEN (20a), 
# & [4 C [2 . 
7, + 01% — = 9=2 N ee 3 5 Fre BEIN RETENE (20b). 


Die Differentialgleichung (20a) für z, ist aber identisch mit der Gleichung (7) für die Staupunkts- 
strömung, da o, identisch mit der Stromfunktion o der Staupunktsströmung ist. Die höheren 
Funktionen r, und z, können daraus schrittweise nach den früher erläuterten numerischen 
Verfahren errechnet werden. 


Der schwingende Zylinder: _ 
Für den schwingenden Zylinder machen wir den Ansatz: 
v—= V, [eos wel + rt? + u + )+sinotig +++). . (21). 
Durch Einsetzen von (17), (18) und (21) in die Grenzschichtgleichung (2) erhält man folgende 
Gleichungssysteme für die 7;;: 


w 
v2 D re 
Tg t AT — a Too = 0 


EEE SE (22) 
nr , o 
Te t 9,9 + - Ta 0 
16 
z 7 (43) ’ GC ’ 
71 + 9791 — A 79 — 2,71 =—3 ., 9 en 


; _ w 7 = [2 . f2 
1. —_— = 
Tat Ta + a 217 20,19 = — 3 a OzTog 


und entsprechende weitere Gleichungssysteme. 


2) Mit den Bezeichnungen von Hiemenz it „= X, 9 =4X, o = 3030 9; für oc, = — 0,597. 
2% 12 


u 

. 
} 
Pr, 
| 
Br 
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Br; ‘ 
j i i i j j i ür die Schwingungsbewe- 
Die Gleichungen (22) sind aber identisch mit den Gleichungen (8) für die Sc) | 
gung bei der ae Das in Bild 3 dargestellte Rechenbeispiel für 7,,; und 7y; ent- 
spricht einer verhältnismäßig langsamen Schwingung mit der Schwingungszahl 


0 BO 

2n An S 
Auch hier können durch Überlagerung von Lösungen mit exponentiellen oder periodischen 
Zeitfunktionen allgemeinere Fälle gewonnen werden. 


N 


6. Zusammenfassung 


Die Grenzschichtströmung an symmetrisch angeblasenen keilförmigen und zylindrischen 
Körpern wird für den Fall untersucht, daß diese Körper gleichzeitig in Richtung der Schnittkante 
des Keiles oder der Achse des Zylinders nichtstationäre Bewegungen ausführen. Ein besonders 
einfacher Grenzfall ergibt sich bei der Staupunktsströmung, wo nicht nur die u-Komponente, 
sondern auch die v-Komponente unabhängig vom Ort ist. Von einem gewissen technischen 
Interesse dürfte die Grenzschichtströmung an einem angeblasenen periodisch in axialer Richtung 
bewegten Kreiszylinder sein. Die Grenzschichtströmung ist als Reihenentwicklung nach der 


Abwicklung des Zylinders darzustellen. Das erste Glied dieser Reihenentwicklung entspricht 


der Staupunktsströmung, die zahlenmäßig für ein Beispiel berechnet wird. 
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The Stress Distribution in a Shaft Press-fitted with a Collar 
. by H.Okubo in Sendai, Japan 


Im folgenden werden frühere Untersuchungen über das zweidimensionale Problem des Pfeilers mil 
rechteckigem Querschnitt auf elastischer Unterlage fortgesetzt‘). Es wird hier das Problem der Welle, die 
durch einen Ring oder ein Rad gepreßt wird, behandelt. Lange hat man von der theoretischen Behandlung 
dieses Problemes trotz seiner technischen Wichtigkeit abgesehen. Die exakte Lösung macht beträchtlich 
mathematischen Schwierigkeiten. Man kann aber auf dem üblichen Wege eine Näherungslösung erhalten, 
wenn aüch die Rechnung etwas kompliziert ist. 


This is the continuation of my previous study on the two-dimensional problem of an elastic foundation 
compressed by a pillar with a rectangular cross section!). In this paper, we treat the problem of a shaft, 
press-füted with a collar or.a wheel. T’'he problem has long been left out of theoretical investigation, regardless 
of its importance from technical point of view. Considerable mathematical difficulties will be encountered 
in finding the exact solution of the problem, but an approximate solution can be obtained by a usual way, 
though the analysis becomes somewhat complicated. 


Dans Varticle suivant, mes &ludes anlörieures sur le probl&me de deux dimensions d’un fondement 
Elaslique comprime par un pillier de profil rectangulaire sont continuses. Nous traitons le problöme de 
l’arbre presse par un anneau ou une rous. Malgr& son imporlance technique, les spe£cialistes ont Evite 
jusqu’4 maintenant & Etudier la th£orie de ce problöme. Pour alleindre une solution ezacte, il faut surmonler 


beaucoup de difficult&s math&matiques. Cependant, on peut oblenir une solution approzimalive par la voie 


usuelle, quoique le caleul soil un peu complique. 


” 


Humecnenymınde ABIAeTcA NPONOMKeHHEM MOHX IpeAbINyINHX HecHe71oBaHHÄ BYXMepHoä 


IPoÖ.IeMsI YIIPyTOTO OCHOBAHHS, CHKATOTO CTONDOM C IIPAMOYTONIBHEIM HOIepe4YHBIM CeYeHHeM. ._ 


B HacrToamei cTaTbe P&accMATpHBaeTcH 3a1a4a BA1a, OÖSKHMAEMOTO KONBIIOM HIIH KOJIECOM. 
JaNaya 9TA, HECMOTPA HA CBO& 3HAYCHHE C TEXHHYeCKON TOYEH SpeHuA, MOATOe BpeMA He 
UOABepFauach TEOPeTHYecKOMy MHocie1oBaHuw. BälABreHHe TOYHOrO PemteHng Ipod.sreM&I 
COUpSIKeHO C BHANHTEJIBHBIMH  MATEMATHYECKHMH TPYAHOCTAMH; IIPHÖNKIKEHHOE pemeHne 
OAHAKO, MO3KeT ÖBITB NOCTHTHYTO OÖBIYHBIM IIYTeM, XOTA AHaJIH3 HeCKOJIbBKO CIIOIKEH. 5 


5 H.Okubo: Quart. J. Mech. und Appl. Math. (1951). 
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Symbols used: 


TIERE Inner and outer radii of a collar. 

b2.2......=0(.Half.length-of a shaft. 

A Stress function. 

0, 08, 0, ... Normal stresses on planes perpendicular to r, © and z directions. 

Vs une Shearing stress parallel to 2-axis on the plane perpendicular to +. 
NR Young’s modulus. 

VE are Poisson’s ratio. 

IRRE: Components of displacement parallel to r und z axes, 

N a RATE Modified Bessel functions. 

ET Tight allowance. 

DER en The pressure on the contact surface of fitted members. 

SEE Mean contact pressure. 

En a en Correcting factor of mean contact pressure for the usual formula. 


Let us consider a shaft of diameter 2a, with its length 2], pressfitted with a collar, as 
shown in fig. 1. The inner and outer diameters and the length of the collar are denoted by 2a, 
2b and 2, respectively. 

We shall take the axis of z along the axis of the cylinder, and r, © as the polar coordinates 
. in z, y plane. When a solid of revolution is strained symmetrically by forces applied at its sur- 


face, the stresses in the body are generally 
) f 
given‘ hy?) 
P) 02x S N 
2 =, 75, RINE 
P) 1 a 2b 
FRAR Eee Z 
Zap >| 2, rör)’ a) 
7) e ’ 
u en „a3 Bar EN 
x Le Rz] 
0) 2a 2 
en ge EN ER 
Geis ja v)V?x FPE) 5 
and the displacements in the directions r Fig.1. 
and z are given in the forms 
a 
NER ERTETT SR: a ar: a2 | r Or 2). 
rovided with | : 
ß A er De TER, a EEE 
where E and » denote the Young’s modulus and the Poisson’s ratio respectively. 
We shall proceed with our calculations assuming that a friction does not act on the contact 
u) surface of the fitted members, and that the normal traction o, which acts on the side surfaces 


does not affect so much upon the pressure distribution on the contact surface ®). Accordingly, 
the boundary condition at the side surfaces of the fitted members can be simplified as 


(£) The shearing stress ,, vanishes. iR Et 
(68) The resultant of the normal stress o,on each side surface vanishes. 


(a) Stresses in the collar 
| According to the above assumptions, the boundary conditions for the collar become 


el; o,=P,+ 3 P„eosnnz, atr=a, 
n=1 


Id, N, 2 atr=bd, a er (4), 


7, =0, 2rjo;rd=0; i ER atz=-+1 


‘where P, is an unknown constant. | 
TIEFER 2 < 6 r A E " . z y 976. h 
2) A/E>H. : The Mathematical Theory of Elastieity, 4th edition (1927) p. 2 : 
s Me en of the latter assumption has been confirmed in a particular case treated in my prece 
ding paper (loc. cit.). Be w 
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If we denote the stress function for the collar by %, then we have 


ed B,2r?+0,2logr+ 5; {A, I, (naar) + BarI,(nar) + CaK,(nrr) 
+ D,rK,(nrr)} sin nz 


where I and K denote modified Bessel functions of the first and second kind respectively, and 
A, B, C, and D, are arbitrary constants. 
By a slight calculation, we have 


(5); 


7, =642+4B,24+2n 5 n{B,I,(nar) — D,K, (nzr)} sinnnz. 
n=1 
Substituting the expression of x, given by (5) into (1), we have the corresponding stresses 
o, and 7,, as 


,=6vh—2 (2) + 
5 {A„Ppın(r) + BuPen(r) + Onden(r) + DuPan(r)} cos nz, (6), 


%, = > {Ankın(r)-+Bykan(r) + Onkan(r) + Dukın(r)} sin nz 
where | 


; le REN araan)), 


Pan (r) = — m? n?{(1—2») I, (nr) + narl,(nar)}, 


Wa 


2 


Pan(r) = — een) + nnK,(nar)| ! 


Pan (r) = n?n®{(1—2v)K,(nar)—nnrK,(nar)}, 
kult)=n’n®l(nzr), 
kon(r) = n?n®{2 (1—v) I, (nar) +narlı(nar)}, 
kan(r)= — nm’? K,(ner), 
 Kan(r) = nn? {2 (1—v) K,(nar) —nark,lnar)}. 
To satisfy the boundary condition (4), we have the following systems of simultaneous 


equations: 
3(1—r) u +2(2—r) = 0, 
6vA—2(1—2r) B, +5 En 


en 1 35 a 
a _ 


A212) +2 =0, 


Pın(a) An + Pan(a) B + Psn(a)Cn + Pın(a)D, = P,, 
Pın(b) An + Pen(b) B. + Psn(b)Cn + Pın(b)D, = 0, 
kın(a) A, + kon(a)B, + kan(a)C, + kın(a)D, = 0, 
kın(b) An+ kan(b) B, + kan(b)C, + kın(b)D, = 0, 
n=1,2,3,... 
Solving the simultaneous equations (7), w. e x 
represented by terms with sr Dh Er A 


(b) Stresses in the shaft 
The boundary conditions for the shaft are 


=), = P,+ 3 P,cosnnz, Maas, 
n=1 ar=a, i 5 
FRbRR 122>1,—1>22—1 SE 
BT en 70. N Or ; 
: an [o,rdr=0, atz=+lI 
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which are transformed into 


Lad, = In 00s"77, air = 
TE NE ee er (9), 
0, 2mfo,rdr—=d, atz=-+] 
where ; 
Po > 
g9=—, NR ER Pe EN 
% m=0 
and 
1 2n(—1)rr! in 9 
Gm=7 when Im=n, Gm Eee we) N, when Im #£n. 
We shall put the stress function x, for the shaft in the form 
g=421Ber ty Anl, “ BT. "sin 2.24.00), 
n=1 


where A, and B,„ are arbitrary constants. 
Then, we have 


V7?,=64A2+4Bor+ ni > BanI, | sin er 
> 2 n=1 
Substituting the expression of %, given by (10) into (1), we have the stresses in the shaft as 


64, —2(1-2)B,4 > {Antınlr) + Bagen{r)} cos "77, 


NTTZ 


o—=6,/h—2(1 2) B+ 2 {An hın(r) + Bahen (r)} a 


0, = 6(1—»v) A, + 4(2—r)B, +2 {Antınr) + Batan(r)} cos nn ö 


2 Bi {Ansın(r) + Basen (r)} sin & E = 


le) 
ne ne rt [ee], 
kn) = — 2). 
wo ul) | 

= Bea) 
ln "a +7 ei, 


- To satisfy the boundary he (9), we have the relations between the unknown con- 
stants An. and By as . x 
2 34 —(-»)B=7; 


ee (12). 
Jin (a) An + gn(a)Bn = 2 
81n(@) I. . — $2n (a) B, ==) 

= 1.238023 


where 


ton (r) = = 
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Solving the simultaneous equations (12), unknown constants A}, B, are represented. by 
terms with coefficient c,, and subsequently by linear functions of unknown constants P,. 


(c) Displacements on the contact surface 


If we denote the radial displacements for the collar and the shaft by,, and u, respectively, 
from equation (1) we have 


ee (2 Bra =.) 11? Si pen2{AyI, (nr) + Bırla(ner) 
n=1 


E E 
— (,K,(nar) — D,rK,(nrr)}cosnnz, EHEN 
2(1+4v) ,. 147 S a ; ka ; nnr\\‘ nnz 
Uy—= ———— Bir — Anl, \——})+ Bar Il, \——); cos 
r E x Eu % ERBE I ı 
On the contact surface, the radial displacements must satisfy the relation 
ö er 
%—urz—0 2 ee ee (4): 


where 6 is the tight allowance, i.e., the excess of the original shaft diameter over the inner 
diameter of the collar. . * 
2 


l 


At the contact surface, cos can be expanded in a Fourier’s series as 


NTz l 5 
cos Tg Mu > CnmCOSMmT2. 
m=1 


Accordingly, from the condition (14) we have a system’ of simultaneous equations for P,, in the 
form 


C : < mn?( , „ [mna 2 Eö 
Br Yanı"zz (Anl("7) + ma2 ("7 ) tag 0. 


A,I(nza) + B.al,(nra) — 0, K,(nna) — D,aK,(n ra) 


1x z MITA x Ma 
az "mn (Anl (®7)+ Bnan (77) =0. 
n=1,2,3, ... 


which can be simplified by virtue of equations (7) and (12), as 
Os. j < mna\., .„ #Bö 
2aB,+ ee 2aB, + en, mond (2) Em Si») =(, 


Be te 


nn 
I(nzo) B, + Rı(naa) D.——, I Mon ("7m 0, 


-. wlan 
The simultaneous equations (15) are-the linear equations of P, and so we can find the approxi- 


mate value of P, by taking first several terms in each series. 
In the particular case where the lengths of the shaft and the collar are the same, putting 


l= 1 in equations (7), (12) and (15), we see that the unknown coefficients P, and subsequently 


Ay, By, On, D,, An and B„ all vanish except the case when n—0, and hence we have the usual 
formula for a thick-walled cylinder as | 


Ps n— er) (b? —.a?) Bun ui Te 2: (16): 


= 


, As an example, we shall consider a shaft of. diameter 2 inches with its length 6. inches, 
which is press-fitted with a collar whose inner and outer diameters and the length are 2, 4 and 
2 inches respectively, then | 


a=1, 'd=2, -I=3, 


1 BERN 
Assuming een from equation (7), we find the unknown constants A,, B,, CO, and D,, 


represented by terms with unknown coefficient P, (The first nine terms in each series are consi- 
dered). From equation (12), we can find A, and B,, represented by linear functions of P,. 
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Taking the first 25 terms, and substituting these expressions of unknown constants into (15) 
we. have a system of simultaneous equations for the calculation of nine unknown constants P,, as 
1.9633 P, +0.0367 P, —0.0187 P, +-0.0107 P, —0.0070 P, 40.0051 P, 
—0.0040 P, +0.0033 P, —0.0034 P, = — 0.3750 E86, 
—0.1752 P, —1.7136 P, —0.0611 P, 0.0358 P, —0.0242 P, 10.0178 1 
—0.0145 P, +0.0126 P, —0.0137 P, = 0, 
0.4481 P, —0.3055 P, —4.7677 P, —0.1030 P, +0.0748 P, —0.0581 P, 
-+0.0506 P, —0.0470 P, +0.0568 P, = 0, 
‚ —0.1709 P, +0.1192 P, —0.0686 P, —2.1774 P, —0.0364 P, +0.0294 P, 
—0.0273 P; +0.0266 P,_—0.0351 P, —= 0, 
0.0841 P, —0.0604 P, -+0.0374 P, —0.0274 P, —1.2366 P, —0.0179 P, 
0.0181 P, —0.0184 P, +0.0261 P, = 0, 
—0.4817 P, +0.3563 P, —0.2317 P, +0.1762 P, —0.1427 P, —7.9674 P, 
—0.1463 P, +0.1485 P, —0.2237 P, = 0, 
0.3153 P, —0.2426 P, +0.1686 P, —0.1384 P, +#0.1204 P, —0.1240 P, 
— 5.5163 P,g —0.1247 P, +0.2112 P, = 0, 
—0.2283 P, 40.1793 P, —0.1337 P, +0.1126 P, —0.1053 P, 40.1056 P, 
—0.1069 P, —4.0109 P, —0.2465 P, = 0, 
0.2025 P, —0.1718 P, +0.1425 P, —0.1317 P, 40.1303 P, —0.1399 P, 
—+0.1602 P; —0.2163 P, —2.6068 P,; = 0. 
Solving the simultaneous equations, we finally find the values of P, as (unit Z6): 
Br =—01913 P= 0.0215 ,=—0.0210, P,=" 0.0183, P,= —0.0161, 
I ° VWOSTLB= —0.0149,-P, =; 0.0157, P, = — 0.0221. 
The corresponding values of A, and B, are given in Table 1. 


Table 1. Values of A,, B,, 


n A, B,, n A, B,, 
0 | 0.2664 - 10-1 0.1598 - 10-1 13 0.9391 - 10-° | 0.8260 . 10-° 
1 2.1119. 10-2 | 0.6150 - 10-1 14 0.1627 .10-° | 0.1444 . 10 
2 0.9001 - 102 | —0.4413 - 102 15 | 0.3559 . 10-10 0.3183 - 10-0 
3 | 0.2401 . 10-2 0.1460 . 10-® 16 | 0.2612 - 10-20 0.2353 . 10-10 
4 | 0.2087 . 10-? 0.1422.. 10-® 17 | 0.4441 . 10-1 0.4024 . 10-11 
5 | 0.2947 . 10-4 0.2155 - 10-4 18 0.1152.. 10-12 | 0.1049 . 10-4 
6 0.2941 .10-5 | —0.2257 . 10-3 19 0.7888 - 10-2 | —0.7219 ..10-12 
7 0.2211 .10-5 | —0.1757.10-°- | 20 0.1288 . 10-22 | 0.1183 . 10-12 
8 0.3611 .10-° | —0.2948.10° | 21 | —0.4075 . 10-2 0.3760 - 10-13 
9 | —0.5446 - 10-7 0.454110” | 22 | 0.2636 . 10-13 0.2441 . 10-22 
10 | 0.4034 - 107 0.3422.10” | 23 | —0.3775..10-14 0.3507 . 10-14 
11 |. —0.8945 . 10-® 0.597510 | 24 0.1997 . 10-1° | 0.1861 . 10-14 
12 0.1271.10-° | —0.1107 . 10-® 


*  Substituting the numerical values of A, and B, into (11), we can evaluate the values of 
the stresses in a shaft. ‚But the series in the expressions in equation (11) do not.converge rapidly 
when r approaches q. The terms taken in each series are still insuffieient for the calculation of 
the contact. pressure, and so the distribution calculated by the equation 


n=8 
2» (0,)r25 = P,+ 5) Nah, COS TUTE Or ee hen Fk Lone Be: (12); 
n=1 a 1 


becomes as is shown by a dotted line in fig.2. For an accurate calculation of the contact pres- 
sure, we must take far more terms in each series, but it makes the numerical caleulations con- 
siderably heavy, so we shall estimate the pressure distribution by assuming that p can be repre- 
sented approximately by the equation®) 


% BelAamWen EN... Zen en (18). 
We shall put a, 
n=8 
erde DEI GEREROENIIT ENTE (19), 
.Nn=l 


4 The appropriateness of thisestimation depends on the suitability of the assumption (18). While the 
en =+ * have some experimental backing. See S. Timoshenko: Strength of Materials, Part II, 
2nd edition. New York 1941, p. 242. Be Sur, i 
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and shall determine the unknown constants @, @’ and k so as to satisfy the condition that the 
integral of e? taken over the whole area of the contact surface becomes minimum, ViZz., 


1 1 2 1 
0 EL: 2 -2|e Br 
lea ee ee ur 
[) 


From which we have 
k 
[4 — 


Q-P+70=0, 


e&+1,, _ k 'S Pfe—1)r—1} 
9 Ar Teer ner aa 
2kek—ek 1 Se NO 


9 - Pot get)? 
k2 SI, Iket(— 1)" A 


 kk—e+1 2 Pu) kKe+n?n? (k®+n? 72)? 


Fig. 2. The pressure distribution on the contact surface. (1) Values, calculated 
by (17), (2) Values, estimated_by (18) 


To satisfy the condition (20), k must satisfy the relation 
2E__02% EN Ph k(—_ 1Ww_ 
(kr—ea+1)|2 ie e?!11 e 55 P„[e*(—1)"—1] 


4(kt—k+1) Ka Ptmm 
Sa 
Be ( N; en. (k?—.n? 7?) al | 1) 
ENT ee 
When e* is very large compared with unity, the equation (21) can be satisfied approximately, if 
"SS P,(—-1)r "SP (—1)r (—n? m?) | 
ES tn 2 BD, (MB +tnn nt ae he Fe (22), 


from which % is determined as 
k—=1637, 

and from the relation (20), we find 
9=—017%6 and = —0.9725 - 10%. 


The pressure distribution on the contact surface, estimated by (18) is also shown in fig, 2. 
From the figure we that see that the contact pressure near the edges becomes considerably large. 
According to the above calculations, the mean contact pressure of the case is 


P,=—0.1918 6, 

while the corresponding value, calculated by the usual formula for a thick-walled ceylinder is 
P,= — 0.1875 EÖ, 

which is about 2.2 per cent. smaller compared with the former, 
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From equation (11), the stresses on the section 20 and also on the axis of the shaft have 
been calculated, and they are shown in figs. 3 and 45). From the figures we see that the predo- 
minating stresses in the shaft are oo and 
o,, and they become considerably large 
at both edges of the contact surface. 
Accordingly, for the discussions of press- 
fit stresses, it becomes essential to find 
the pressure distribution on the contact 
surface. Moreover, if only it was found 
the further analysis could be worked out 
by a known method. 

It becomes sometimes of technical 
importance to find the substantial mag- 
nitude of the mean contact pressure. 
In practical calculations, however, the 
end effect is usually neglected and the Fig. 3. The stresses on the section z = 0 
formula for a thick-walled cylinder is 
used. But we shall investigate the problem by taking the end effect into consideration, and 
a correcting factor for the usual formula will be derived. 


Stresses 


As is seen in the preceding example, the coefficients P,, P,, ..., are smaller quantities 
compared with P,, so neglecting them, we have (approximately) 
mrca Gmo Po 
: | re RN 
Bat; j ) m? I F a ee 23) 
I ae er mar | 


Stresses 


Fig. 4. The stresses on the axis of the shaft 


Substituting the above relation into the first equation of (15), we have 


a(b?--ka®) ee ERENDME SE 24 
er I an“ Bauen Ei 
where s 
. sin? 

en m —_ 5 re 

Sk BEE TE ERENT E re e e  W and k=——, 

na nna\ 2%2(1—») Em 

Zu ne een 


from which the approximate value of P, can readily be obtained. 
The equation (24) coincides with the usual formula in the particular case of I=1, and for 
a shaft ei vokar : the dimensions in the preceding example, P, becomes —0.1910 Eö which 
almost agrees with the previous result. 
Putting ?, in the form 
EBö (b?—a}) 
ne un (25), 
4ab? 


the coefficient & for various dimensions of a shaft and collar has been calculated from equations 
(24) and (25). The results show that & increases according as l increases, but remains almost 


I 


5) The value of og:on the contact surface is obtained from the relation =, +77, 
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constant when I becomes larger than 1.5. Accordingly, the values of & only for I=a are shown 
in Table 2, » being assumed to be 0,27. 202 

In practice we often meet with the case where a shaft and collar is of different material, 

e. 8:, the former is of mild steel and the latter is of cast iron. The analysis is, however, quite 

similar to the case of same material, and instead of 


equation (24) we have 


Table 2. Values of & 


au ah kz 4(l—m)hUS:|n 
a=1,1=2 | 1.02 | 1.02 | 1.03 | 1.03 2 a2 l anı 2 
a=2,1=2| 1.02 | 1.04 | 1.05 | 1.07 3 (26), 
a=3,1=3 | 1.03 | 1.04 | 1.07 | 1.11 Be RE _y 

a —=.4, i=4 1.03 1.05 1.07 1.14 2 (1-49) > 

a=5,i1=5 | 1.03 | 1.05 | 1.08 | 1.17 


where E,, E,, v, and », denote the Young’s moduli and the Poissons ratios of the collar and 
the shaft respectively, and 


ee lt E,(1+9)' 
nn EEE 
x ner 
SA > . 
— NTA NTTa 212 (1—»,) 
le Yr( Eee 
Thus, the correcting factor for the case is given in the form 
€ Ed 
PR=—-—n nn RUHE. 
Fa: EN 
20, A E, 


Assuming the elastic constants of cast iron and steel as Z,—=8000 tons/in?, E,—= 13000 
tons/in? and »,—=9,—=0.27, the values of &’ have been calculated from equations (26) and (27). 

According to our calculations, the values of the correcting factor for the latter case are 
somewhat smaller compared with those of the former case (the table for &’ is omitted here), 
and it shows that the value of a correcting factor becomes smaller with the increase of the ratio 
of the Young’s moduli of the shaft and the collar. 
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Kugel mit normalgerichteten Einzelkräften 
Von ©. Weber in Pullach bei München 


‚ Untersucht wird die Kugel mit Einzelkraft normal zur Oberfläche. Die auftretende Belastungsfunktion 
wird durch Inversion bestimmt. Durch Überlagerung findet man den Belastungszustand durch mehrere 
Kräfte. Näher untersucht wird die Kugel mit zwei an den Polen wirkenden Druckkräften. 


64 single force is acling upon a sphere, normally to the surface, and the stress function is determined by 
inversion. If several forces are acting, Ihe state of stress is found by superposition. The case of two forces 
at the poles of the sphere is investigated particularly. 


‚L’auteur examine la boule & force single agissant normalement & la surface. La fonction de chargement 
Er se prösente est döterminde par inversion. Par superposition on trouve l’tat de chargement par plusieurs 
orces. La boule avec deux forces de pression operant aux pöles est soumise & une recherche en detail. 


- Hecnenyerca chepa ce onumounoli cmıoli, HOPMaıbHoH IOBEPXHOCTU. Bostmkamınaa Pyakuus 
HATPy3KH OTPeNenaeTchH TOCPeACTBOM ea Ilyrem nepekpsrrus MokeT ÖBITb HalieHo 
COCTOAHH® HATPY3EH HECKOABKHMH cHıaMHu. [lonpoöHee Hecnenyerca chepa c AByMA CiKHMaI- 
INMHMH CHIAMH, MEHCTBPIOIMUMM Y_ IIOIOCOB, FAN Fr 


„Mit Inversionen und anderen Kniffen“ 
1. Einleitung 


„ Wirkt auf eine Kugel aus isotropem H o o ke schem Stoffe eine äußere Belastung, so gibt 
es ein allgemeines Verfahren zur Berechnung der Verschiebungen und Spannungen!). In vorliegen- 


zo 


!) Siehe Love-Timpe, Lehrbuch der Elastizitätslehre, Kap. XI. IE BRLE 


richteten Einzelkräften 2 187 


n are inen rotationssymmetrischen 
Fall vor uns haben. Für diesen Fall wird künstlich. Gleichgewicht hergestellt. Bei der Unter- 


4 suchung wird die auftretende Belastungsfunktion durch Inversionermittelt. Weiter werden - 
Z Reihenzerlegungen rotationssymmetrischer Potentialfunktionen dadurch vereinfacht gefunden, 
E- daß sie erst nur für die Rotationsachse bestimmt werden. Den Fall mit mehreren Kräften. die 

im Gleichgewichte sein müssen, findet man durch Überlagerung. 

__ Als Beispiel wird die Kugel mit zwei Druckkräften an zwei entgegengesetzten Polen unter- ae 
Be ' sucht und die Spannungen im Mittelpunkte berechnet. Die Berechnungen der Spannungen im : 
ER . Randgebiete müßte mit Hilfe schlecht konvergierender Reihen erfolgen. Diese lassen sich ver- a 
Er: bessern, indem man geeignete Belastungsfälle abspaltet, die geschlossen angegeben werden können. Be . 
Es kam mir besonders auf die Anwendung z.T. neuer, z.T. alter Verfahren an. Ich hoffe, daß ve 

nicht allesschon bekannt ist, was man ja nie wissen kann. 4 


air 


* 


ri = 2, Lösung für die Kugel mit rotationssymmetrischer Belastung 


Wir wählen Kugelkoordinaten 7, , ® mit dem Koordinatenanfangspunkte im Kugelmittel- ; 
punkte, bzw. @yz-Koordinaten = rsin d cos 9, y=rsin dsin 9, 2=rcos®. Dann drücken wir die 


PN 


IND 


Zerrungen durch die Ableitungen der Verschiebungen u, v, waus: 


E 
Br: 
Br. 
2 # 
E* 


. 


Er 2 & =; und zyklisch, y,, = . > 


und bezeichnen die Raumdehnung mit &,: 


ren Se oe at: 
R Kacy BE + eu + ..bzw.: a a a (1). 


0: - Nachdem Ho o keschen Gesetz werden die Spannungen N SE NEE 


nenn 


2 a / 2 Zei aner r. Uu.Zz. = - En. 


\ x 
= x = Eu z n 


r diese Ausdrücke für die Spannungen in die Gleichgewichtsgleichungen ein, so T 


_ 


en ER. 2. a ee 
; N A ee, RR Bin 1-2 y: 8 ET rl ” a ee (: ) \ 


ee er NE £ er Er NR ER 
leichungen dieses Gleichungssystemes werden nach &, bzw. y, bzw. z differentiert 


a N en 
a 


2} IE u=pMmTt Pen ir a er Fe tt a . ‘). 
friedigt das Gl.-System (3), 


FF 


ER Se TE s en 


em (5) oline u, 9, Pu stellen Teillösungen 


ist der Ansatz eine vollständige Lösung. Setztmanu,v,w 


x 


en verknüpftsind. 
e Ansätze möglich, bei denen 
Bei der Untersuchung einer 


. ” '. * “ h. 
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Nun kommen wir zu den Randbedingungen: An der Oberfläche der Kugel vom Halbmesser 
»— asind die belastenden Kräfte gegeben. Am Oberflächenteilchen d0 greifen die Kräfte 8, d0, 
H,d0 und Z,dO an. Dann gelten die Gleichgewichtsbedingungen der Oberfläche: 
(805 + Ylay-+ 2Tzz)r-a= 0a UZ. 
Führen wir o,, 7,, usw. nach G].-S. (2)ein, so erhalten wir nach Teilung durch @: 


ou ou ou ou 0% ow 29 _ U. 
rer er 


Weiter setzen wir w, v, wnach Gl.-S. (ö)ein: 


ee. 
\errmeer ae tray n20 (8) 


Oy Oy 
ar bei ER Na ER, ‚zZ. 
+2r 5 +2ar, +2r al ga U:2 


Wir wollen nun zeigen, wie die linke Seite als Randwerte einer Potentialfunktion dargestellt 
wird. Hierbei verwende ich folgende Beziehungen: 
Ist 9, eine Potentialfunktion, so ist auch 


991 091 9 _ „9ı 
N De ee ee 
eine Potentialfunktion und ebenso 
9pı _ ,9Pı 
Oy Fan 


In Gl.(8) stellt die erste runde Klammer die Potentialfunktion „Be dar. Für die zweite 
‚Klammer schreiben wir: 


nee 
(VE Near er FINE EI de 


(v2) 4 (‚me .)_o| ay 
- (y rn + 2 —2/(1 2 )y+ß—4r)r : 
Damit gibt Gl.-S. (8): ö 


Pu | 0, er | Pu = Op Oy 
[ a Ur Ta I a re Br Ar? 20] 


0) &4 a 

2ar.r — — =-E N % 

| inzE "ar dnlr=a Tee 
Auch die eckige Klammer ist eine Potentialfunktion, was durch Differentiation geprüft 
werden muß. Nur das letzte Glied der linken Seite ist in dieser Form keine Potentialfunktion. 
Setzen wir bei diesem Gliede jedoch r? = a®, so nimmt es die Gestalt 20 7 an und wird da- 

| r 0% 
mit zur Potentialfunktion, sodaß auf der linken Seite jetzt nur Potentialfunktionen stehen. 
Weiter formen wir dieses letzte Glied um in 


o[ y n 
a ee! 
2a .(r or v) j 
und erhalten: 
Pu u Op ( 09 ®) 09 09 
en en —u EN EI RE EN 
( eye +2) + 0% a +(e 0% 29) 


r 2 an (9a), (9b), u. (9e). 
ee u ie Te Be 
+| "7 2ay—4ar ar +20? (rt v)=4 a 
Dieses Gl.-S. gilt fürs erste nur für die Kugeloberfläche. Auf der rechten Seite d i 
chungen standen die Oberflächenbelastungen £,, H, und Z,; wir nehmen an, daß sie die Bes 
werte der Potentialfunktionen 8, H, Z sind. Dann gelten die Gl. (9a), (9b) und (9c) für die 
gesamte Kugel. Zur Bestimmung der vier Potentialfunktionen Pus 9% Pu und Y stehen uns 
a die = Gl. (7), Ge). (9b) und (9c) zur Verfügung. 7b AR | 
Nun wird es zweckmäßig sein, aus diesen vier Gleichungen eine Gleichung für nur ei 
Funktion aufzustellen. Wir wählen hierfür die Funktion Y, da sie mit ke Bönalicaene 


ER 
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richtungen verknüpft ist; yist in Gl. (9a) mit & multipliziert, in Gl. (9b) mit y, in Gl. (9c) mit 2. 
Folglich versuchen wir, Gl. (9a) nach x, Gl. (9b) nach y, Gl. (9c) nach z zu differenzieren und die 
drei erhaltenen Gleichungen zu addieren, um eine von den Koordinatenrichtungen unabhängige 
Gleichung zu erhalten. Dieses ist de gedankliche Begründun g für das im übrigen 
bekannte Verfahren. 


Die Glieder mit @,, 9, 9u heben sich hierbei fort bis auf den Ausdrück 


Pu , IP a 
Gere), 


- den wir nach Gl. (7) durch einen Ausdruck von yersetzen. Damit bleiben auf der linken Seite 


nur Glieder nach, die yund die Ableitungen von y enthalten; die Gleichung lautet: 
0% . Op a E oH 92 ) 

che, BA m ee ee Aare ae 
errang, yı% 


Bis hierher hielt ich mich an das bekannte allgemeine Verfahren. Nunmehr will ich 
annehmen, daß eine rotationssymmetrische Belastung vorliegt. Für diesen Fall läßt sich die 
weitere Entwicklung vereinfachen. Z ist eine rotationssymmetrische Potentialfunktion und 


02 m BE : 
ebenfalls — , Zund H sind die Ableitungen einer rotationssymmetrischen Belastungsfunktion R; 


02 
Sr 
0% 0 
N i : r oD Fr OD 5 
Auch » und 9, sind rotationssymmetrisch, während 9, = Fri und 9, = = als Ablei- 
tungen einer rotationssymmetrischen Potentialfunktion ® dargestellt werden können, 
Alle rotationssymmetrischen Funktionen drücken wir entweder durch die Kugelkoordinaten 


r und ® oder durch die zwei Koordinaten z=rcosd und =J® + y=rsin® aus. Die 
Belastungsfunktionen Z und R zerlegen wir weiter in Reihen: 


Zeit Gr) + rel) + *:° = 2 „Cent" 9 (9) 
RB AH) +: Ze rntı Hd). 


Bis auf die Koeffizienten c,, und c,,„»+1 erhält man die Teilfunktionen r"g,(d) auf folgen- 
dem Wege; wir setzen zuerst 
: r9@) =1. 
Durch Integration nach z erhalten wir: 
1 
1l 
Durch nochmalige Integration nach 2 und durch Beachtung, daß das Ergebnis wieder eine 
Potentialfunktion werden muß, erhalten wir: 


Hg(d)=1j2= 1770089. 


r29,(9) = Ie—Irt=, 12 cos W—rsind. 


2! 
In gleicher Weise finden wir die weiteren Funktionen. Allgemein wird: 
Dee IE mean 
NE er Tre r TERN? r BAR 
Die Richtigkeit ist leicht zu prüfen, indem wir den Ausdruck 
en. + t 2 y' ı yn (9) 
oe vater) I 


“bilden, der gleich Null werden muß. 


Die rechte Seite der Gl. (10) gibt: _ 


2 (BE, 0)__ 0 (er er, 22) 
tn BT dat Toy 
YA )=-5 
ala RZ. 


(an — an) md). 


9 


a 
> “ 2 er ne TEE 
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Aus Gl. (10) folgt, daß wir hierzu 
2 "m (d) 


n=0 
setzen müssen. Führt man diesen Ks für yin Gl. (10) ein, so erhält man: 
5 Gp Enke # (if 
Mr ER, 10 Ra — 1%? +2 N) +++) 
era ist y bei gegebenen Z und R gefunden. 


Zur Bestimmung der noch fehlenden Funktionen 9, und ® benutze ich Gl. (7) und Gl. (9c), 
letztere, da sie. symmetrisch in & und yist. Gl. (7) gibt: 
2 2 
rate 2 5 MB Ann Dlmar- 194.10) 
„0% n= 0,1 16% 
mit 
an A 
oa? | ay® 02? 
folgt nach Integration nach z: - 
od 
Sa y=-— ae £ (1-29) +8 —4r)n —1)]n-ırm(d) - - » (12). 
Gl. (9c) lautet: 
„Pu ( 0 ® 02® ( 0? ® =.) | Op u =] 
De EI Fe Dh A, EEE Kr aa 
a NE Baer} De) aa men 
0 


Für die rechte Seite schreiben wir: 


72= PN 9n(9) = Fi ES 
Die Glieder der en Seite, die Ableitungen von ® enthalten, geben zusammen: 
Br a en. u 
0% 02 0y 02 02: ua 0% 02 0y 02 ? dar Or 92: 


Die Glieder in der eckigen Klammer wollen wir getrennt untersuchen. Hierzu nehmen wir den 
Ansatz von y: r' 


ya Hd) EI reanı gulli, ) 


und beachten nur die Glieder mit der höchsten Potenz von 2. Dann wird: 


, 


2 


ICE aan A: )=20n. A al.) 
— 22y = — 2220 a1. -) = RER ar —:- ) | 
SuAzr = —Nlznı re pt a) ran... 


Die einzelnen Glieder sind keine Potentialfunktionen, wohl aber ihre Summe. Darum gibt 
die eckige Klammer den Ausdruck: 


N 2 ET 
— 23 Bale 2: ) = —2 Inn, (9) = —2 N ne nr 79.9)? 
Nun.bleibt noch das letzte Glied der linken Seite nach, dieses wird: 


Or 2 Or 
Setzen wir alle diese Ausdrücke ein, so erhalten wir die- ‚Gleichung: N: 


9y, 008 a 
ar nF PLN +2 Zaun) 


2a nn le): 


(a 00 Ei 
202. (r?_y) 2, 2 Nm NEL PeR RE BE 
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Bei allen Gliedern tritt der Schlußoperator r 2 auf, den wir fortlassen können; wir 
erhalten: or 


od a 1 | 
22 + Eu => [& Er, 1" 9,(0) +2 m-ı Nr 9,(0) — 2 m art 2 9): 


Diese Gleichung gilt für alle Werte n>0. Im Reihenansatze für Z nahmen wir auch das Glied 
6z0; dieses Glied gibt aber eine resultierende Kraft in z-Richtung, da 
1,40 = 6 0 40? 60 18. 


Alle anderen Glieder der Reihenentwicklung geben keine Kraft in z-Richtung. Wir sehen, 
daß bei einer Belastung, die im Gleichgewichte ist, das Glied mit C,, auch nicht auftreten wird. 


k 0 
In der erhaltenen Gleichung für = + 9, können wir im letzten Gliede noch n— 2 durch 
n ersetzen und wir erhalten: 2 


98 1 
rg, = a erahnen la) 
02 @G N 


n=1,2,3... 


Gl. (12) gibt — =  9,.61.(13) = + 04 hieraus folgen = und 9,. Um ® aus = zu bilden, 


_ ersetzen wir nur r*q,(®) durch ‚r+1 „r1(0). Hiermit sind alle Funktionen, die wir für die 
9. In+ 


Lösung benötigen, gefunden. 


3. Kugel mit einer Druckkraft P im Pole z= «a 


Zwar besteht kein Gleichgewicht, das hindert uns aber nicht, die Belastungsfunktion Z 
zu finden. Die Belastungsfunktion R wird gleich Null, da keine Kräfte senkrecht zur z-Achse 
auf die Kugel wirken. 42° 

Um Z zu finden, betrachten wir erst nach Bild 1 i 
den Halbraum 2*>0. Im Koordinatenanfangspunkte ' 
+=y*—=z*=0 wirke die Druckkraft P* nach oben. 
Die Randbelastung sei gleich den Randwerten einer drei- 
dimensionalen Potentialfunktion Z*. Diese muß rotations- 
symmetrisch und von minus zweitem Grade sein, 


damit das Integral Z*2rnr" dr*.über eine beliebige Bild ı 


0 5 RE : 2* 
Fläche 2* — konstant stets denselben Wert gibt. Wir wählen darum hierfür die Funktion c* a: 
Dann wird 


f * 2* 2 g d +’ 2 * [ Bar (7) 5 2 7c* 
c* — anı" ar = 2700” Tue! 33 
Jazer IN ee N" 


Denselben Wert erhalten wir auch für z*— 0; das ist der Wert der gesamten Randbelastung. 


Rape MER, 
Diese gibt die Kraft P*, so daß DEZ wird. ‚Hiermit 
wird 5 


p* 2* 
+ — I Da 
Are In rr3’ 
auf der z-Achse wird 
P* 
Zw = == Duni zrr er a DT er Vase | (14). 


‘Nun legen wir den Koordinatenanfangspunkt nach 


*—— 5; für diese neuen Koordinaten nach Bild 2, die wir ebenfalls mit x*, y*, z* bezeich- 
a > 
"nen wollen, wird 4 
SE: nr 
gi De 
2 i 1 2 0, A 
| ER 
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oder in Kugelkoordinaten 


1 
r* cos #* —— 


P* 2a 


r* 25008 0 + 93 


7* 


1 : e 
Auf der Oberfläche des Halbraumes 2* > 32 ist überall Z* gleich Null, bis auf den Punkt 
” = 0. . . 
\ Nun bilden wir den Halbraum durch Inversion auf die Vollkugel nach Bild 3 ab, 
indem wir r* = 1/r setzen. FG 
In der Funktion Z* setzen wir überall 1/r an Stelle von r* ein; um aber wieder eine Poten- 
tialfunktion Z zu erhalten, müssen wir noch Z* durch r teilen. Wir erhalten hiermit: 


1 
Fo usa _ P*(2a)? 2arcos# — r? 


er 2n [4a?—Aracosd +?” 


Ze na m ee Ta. 
BT E 
(4) Ph er 2a 
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Br 
T 


_P*Ra) Qa—2)2—r” 


2r [2a — 2) +r?] 1? 
Die Kugel, Bild 3, hat den Durchmesser 2a, der Koordinatenanfangspunkt liegt im unteren 
Pole der Kugel. Wir verschieben noch nach Bild 4 den Koordinatenanfangspunkt in den Kugel- 
mittelpunkt, indem wir 2+.« an Stelle von z setzen 
und erhalten: 
Z- 2aP ® — 2? —r'? 
a e=irrpe 
Auf der Oberfläche der Kugel gibt Z überall den 
Wert Null bis auf den oberen Polz=a. Hier wird 
Z= ©. Setzen wir r’=0 und lassen wir z gegen a 
gehen, so geht 


N 
n (a—.)? 
x h | Pr tal 
2 In der Nähe dieses Punktes muß entsprechend Gl. (14) aber Z gegen —— —— — 4_ 
> gehen, um dort die Druckkraft P zu geben. Hiermit wird endlich 2n(a— 2)? 


Be": ( P a? — 2? — r’?2 
E ; Z= SA [(a — 2)? + r?pR 00 SC er ee a ES (15). 


Be’ Bild 3 Bild 4 25 


| Wir sehen, wie durch Inversion die gewünschte Belastungsfunktion gefunden wurde. 
Ber. Es handelt sich hierbei aber nicht um die Inversion einer Spannungsfunktion, also einer Bi- 


: potentialfunktion, sondern um die Inversion der Belastungsfunktion, die ja eine Potential- 
ee - funktion ist. 


Fi Die gefundene Belastungsfunktion Z ist jedoch noch nicht in Form einer Reihenentwicklung 


I, dargestellt. Um diese zu erhalten, setzen wir in Gl. (15) erst r’= 0 und entwickeln Z in Reihe 
Bi. nach z: Rt 


E Du an lit lır33+3(@)+)- 


Be, pP} . 
E. amlirsa tl) +), an+olz]. 
- Damit erhalten wir Z selber, indem wir z* durch n! rg (0) ersetzen: SER 
a er tYa(2) 00) 
mit n= 0,1,2... Die Untersuchun 
zeigten. Nur das Glied mit n = 
- Glied 


Z=— 


8 wird weiter fortgesetzt, wie wir im Abschnitte 2 allgemein 
0 darf nicht auftreten. Wir fügen also unserer Belastung das 


PIE: } AR 5 
1 hinzu. Dieses ist eine gleichmäßig auf der Oberfläche verteilte Belastung in 
Richtung. Auf die Kugel wirkt also jetzt im oberen Pole eine Kraft nach unten und auf der 


u 


- Nach GI. (21) und (22) folgt: 
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Oberfläche eine verteilte, nach oben gerichtete Belastung 


Ange’ so daß Gleichgewicht besteht. 


In der Reihenentwicklung von Z nach Gl. (16) müssen wir dabei n — 1, 2, 3 usw. setzen. 

Die ‘so gebildete Belastung Z kann trotzdem benutzt werden, wenn auf die Kugel die 
Normalkräfte P,, P, usw. wirken, die im Gleichgewichte sind. Für jede dieser Kräfte bilden wir 
die Belastungsfunktion Z; außer den Kräften wirkt dann noch auf das beliebige Oberflächen- 


teilchen dO die Belastung u dO entgegengesetzt P,, en dO entgegengesetzt P, usw. Und 


diese Belastungen sind im Gleichgewichte, d.h. sie heben sich fort, so daß nur die Einzelkräfte 
nachbleiben. 


Wirken auf die Kugel zwei Druckkräfte am unteren und oberen Pole, also in den Punkten 


2=—aund z=--.a, so wird die Belastungsfunktion 
: | a a? — 22 — r'2 
Ser Ana? (a — 2)? + r’27]212 Tr? (+2)? + rapel| "209 (17). 


Die Reihenentwicklung hierfür lautet: 


2=— 8117 9)+7: 312) +) ng tn, 


“Diesen Fall, also die Kugel mit zwei gegenüberliegenden Druckkräften, wollen wir als Beispiel 


weiter untersuchen und hierfür alle Funktionen bestimmen. 


Wir fanden: A | 
Z= a "In emit - ,> SEERTTE QArn-+1)ni..... (19) 
R=0 
r 2% (An —+1)n! 
rer —1 ES 

v3 mar md) mit mid MP Fa Fo) EG») (20). 

Hiermit wird nach Gl. (6): 
eG ar ind) mi, ne al —2)AaRr — 1) 5.2. (21). 


Nach GI. (12) wird: 
od 7 (2n + 1)n! [1 +29) + (3 —4rv)(n — 1)]/r\r 

% ar ee (r— 1%? +2 (rn —1)1+»v)+(1+») & M(9) (22) 

nach Gl. (13): 


(2n-+1)n!n 


Rn) P | Ee BER IEN ERNA AR ar En 
en, | Ant MONFRETFIE-DArNFUH) (23) 
...—.. 0, 2% | 

R +1? +2 HF YA+N)+AFn\a} 


4. Berechnungen der Spannungen im Kugelmittelpunkte 


Grundsätzlich bietet die allgemeine Berechnung der Spannungen keine Schwierig- 
keit. Ich will für die Kugel mit zwei Druckkräften an den Polen die Spannungen im Mittel- 
punkte M der Kugel berechnen. Alle auftretenden Größen versehe ich darum mit dem Index M. 


Da für den Mittelpunkt r = 0 ist, genügt es, alle Größen für » = 1 zu bestimmen. 


h GI. (20) wird 
Nac (20) ä 2 


Ne! + v)na2G 
und hieraus nach Gl. (6) ® Tr Bes : 
3 ua 4(1+9) na?@' 


od st-2ar) P- 
S. = De), = 441+ v) ned” 

B 2. ee ne APe 
Bere stz neue 
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EEE A era ee 
Hieraus 
ar.) 51 3 Se 
(2 v„ U2A+» 2(9+59)2na?@ 
und 


od | 3» 15 | Ba 
(=  Aa-+») 2(9+5»I2naG 
Aus letzter Gl. folgt: 


u) 0?® 1 , = 34 v =E 15 I on P 
(u (2 022 2(1-+%») ' 4(9+ 5») Ira‘ 
Nun finden wir die Spannungen, wobei stets im letzten Ausdrucke » = 0,3 gesetzt ist. Aus £; 
finden wir a 
Ei 3-P 
Br - ra 7 3» gi Ehe: 
; Weiter ist: 
| sale a Be 
ee en - 2 2(9+15»)/na? we T 
und 


v 15 5 we: 
es a a Je 
v-20(; er 29 e) 4(9 + 15») wa? ed 
Im Mittelpunkte treten folglich in z-Richtung Druckspannungen und senkrecht dazu Zug- 
spannungen auf. 


5. Weitere Berechnung der Spannungen 


Wollen wir die Spannungen im Randgebiete berechnen, so konvergieren die Baer nur 
langsam. Die Pole sind Randpunkte des Konvergenzgebietes, wie Bi weitere Begrenzung dieses 
Gebietes aussieht, werde hier nicht weiter untersucht. 

In den Polen selbst verhalten sich alle Größen wie bei einem Halbraum, der durch die 
Druckkraft belastet ist. Für den Halbraum mit einer Druckkraft P im Koordinatenanfangs- 
punkte lautet die bekannte Lösung: 4 


3P2 Pia | P( 
a a en 
REP IL 2) Tu 3 2 a 
; MEERE en re) 


BEINE en 1 e -) 
9a TI N 
Um nun den Zustand der Kugel in Randgebieten zu untersuchen, 
zerlegen wir ihn in die Zustände, hervorgerufen durch die Druck- 
kraft P im Halbraume z<.a, durch die Druckkraft im Halbraume 
12x “ 2>—.a und den Restzustand; Bild5 veranschauliche die zwei 
H \ “ersten Teilbelastungen. Wir wollen nur (0,4 0,4 0,) weiter ver- 
u ne 


folgen. 
Für die zwei ersten Anteile erhalten wir aus Gl. (24): 


Be N) a+2 a—2 
NE Imrertme 5* | 


Nun müssen wir diesen Ausdruck in Reihe zerlegen. Wir benutzen 
wieder die Möglichkeit, die Zerlegung für. = V9 durchzuführen: 


let pr Ho | El or | 
@R ra ee )=- DP(+n) | ale) ” 


ar BE 


Bild 5 


n=1.3. 


Folglich wird allgemein: | Bern | | 
SE ; 2Pp(: TER EV 
177 als 0, ,—= are 3 n! (2) äh EEE, $ (25). 


“ 
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Unsere Lösung für die Kugel gab nach Gl. (20) und (21): 


%+r%,+,—=2G8 Bern 
1— 2» 


AP») [n + 3) 


5) 
2 


n® „zn —1)+2(r—1)(1+v)-+(1-+») a1) Mm—1(d) (26), 


Für große Werte » geht der Bruch hinter dem Summenzeichen gegen eins u | 
nd I 
der Ausdruck der Gl. (25) gegen den Ausdruck der Gl. (26) Eh: : re 


Die Lösung für o,+0,-+- 0, der belasteten Kugel können wir schließlich in folgender 


Form schreiben: 


sa? 


te] a+2 +7 
r 


bee 


ar 
(27). 


2) 


re rn - m) no) 


n=1,3... 


Die jetzt auftretende Reihe konvergiert offensichtlich besser. 


* Eingegangen: 7. Dezember 1951. 


BUCHBESPRECHUNGEN 


F. Strecker, Praktische Stabilitäts- 
prüfungmittels Ortskurvenundnu- 
merischer Verfahren. XI-+1898. mit 
101 Abb. Berlin-Göttingen-Heidelberg 1950. Springer- 
Verlag. Preis brosch. 15,— DM, geb. 18,— DM. 


In einer Monographie ‚‚Die elektrische Selbst- 
erregung mit einer Theorie der aktiven Netzwerke‘, 


—(Ref. diese Z. Bd. 28 (1948) S.283) hat Verf. die 


Stabilitätsuntersuchung schwingungsfähiger Systeme 
nach der Ortskurventheorie ausführlich begründet 
und dargestellt. Das vorliegende Buch soll diese 
Methoden einem weiteren Interessentenkreis zugäng- 
lich machen und dem praktisch arbeitenden Ingenieur 
eine einfache Arbeitsanweisung zur Behandlung auch 
schwieriger Stabilitätsfragen liefern. 

Die grundlegenden Betrachtungen sind, weitgehend 
auf die elektrische Verstärkertechnik zugeschnitten, 
fast alle Untersuchungen knüpfen an die Frequenz- 
gänge geläufiger Übersetzungsfaktoren, Scheinwider- 
stände usw. der vorgelegten Systeme an. Zum Ver- 
ständnis der Vorgänge ist jedoch eine Erweiterung 
der in der Wechselstromtechnik üblichen Methoden 


. nötig, die aber insbesondere durch geeignete Be- 


zeichnungswahl der gewohnten Denkweise des Elektro- 
ingenieurs angepaßt wird. Dabei werden physikalische 
Vorstellungen bevorzugt; auf eine vollständige 
mathematische Durchführung wird zugunsten von 
Gebrauchsanweisungen, Rechenschemata und zabl- 
reicher Beispiele, besonders aus der Hochfrequenz- 
technik verzichtet. Statt mathematischer Vorkennt- 
nisse wird physikalisches Einfühlungsvermögen ver- 
langt. 

Er Stichworte zur Kennzeichnung des Inhalts: 
Grundlagen der Stabilitätsprüfung mittels Ortskurven, 
Begriffe und Bezeichnungen der Netzwerkstheorie, 
Umlaufs- und Fortsetzungkriterien mit numerischen 
und graphischen Verfahren und zahlreichen durch- 
gerechneten Beispielen, Kriterien für den Stabilitäts- 
grad, nichtlineare Systeme (insbes. Pegelwandler). 

“ Mitdiesem Buch wird in derumfangreichen Literatur 
über Regelungssysteme eine Lücke bezügl. der Orts- 


kurvenverfahren geschlossen. — Es fehlt allerdings 
“auch hier (neben einigen neueren Verfahren) eine Dar- 


stellung des in der Begründung anspruchsvolleren in 
Amerika entwickelten ‚‚Bode-Verfahrens‘“. 
Dresden.’ -- N.J. Lehmann. 


R 


Jenaer Jahrbuch 1951. Wissenschaftliche Veröffent- 
lichungen des Zeisswerkes. Herausgeber: Optik Carl 
Zeiss Jena VEB. XII -+ 2908. m. 174 Abb. Jena 
1951. Kommissionsverlag Gustav Fischer. Preis geb, 
20,— DM. 

Geh.-Rat Harting, der die Herausgabe dieses Jahr- 
buches betreute, ist vor seinem Erscheinen gestorben; 
so eröffnet den Band ein Nachruf auf Harting. Im 
übrigen enthält er 20 Abhandlungen, die teils histo- 
risch orientiert, in der Hauptsache aber Problemen 
gewidmet sind, die sich im Forschungsbetrieb der 
Zeisswerke ergeben haben und die fast alle von 
wissenschaftlichen Mitarbeitern dieser Werke bearbei- 
tet wurden. Von diesen Aufsätzen werden vielleicht 
besonders die folgenden drei die Leser dieser Zeit- 
schrift interessieren. Focke, Realisierung der Herz 
bergerschen Theorie der Fehler fünfter Ordnung in 
rotationssymmetrischen Systemen mit einer Anwen- 
dung auf das Schmidtsche Spiegelteleskop, die an die 
Arbeiten von Carathöodory über geometrische 
Optik anknüpfend mit den von ihm entwickelten 
Mitteln Probleme der praktischen Optik angreift und 
ein Verfahren entwickelt, das bei Anwendung auf die 
Berechnung des Schmidtschen Spiegelteleskopes im 
wesentlichen zu den gleichen Ergebnissen führt wie 
eine von Carath&odory angewendete spezielle Methode. 
Kneschke, Die zeitlich veränderliche Abkühlung 
von ebenen Kreisplatten und Bauer, Wärmespan- 
nungen in konischen Scheiben. — Möge es möglich 
sein, daß diese Jahrbücher, dessen erstes 1950 er- 
schien, weiterhin regelmäßig von den Zeisswerken 
herausgegeben werden. 


Dresden. i Willers. 


Prof. A.Rohrberg, Theorie und Praxis 
des logarithmischen Rechenstabes 
(Mathematisch-Physikalische Bibliothek, Reibe I 
Nr. 23). 9. Aufl., 648. m. 15Abb. Leipzig 1951. 
B.G. Teubner Verlagsgesellschaft. Preis 1,40 DM. 

Die in der Z. angew. Math. Mech. 30 (1950) 8. 334 
angezeigte achte Auflage wurde so ergänzt, daß jetzt 
die in Deutschland heute gebrauchten Rechenstäbe 
„System Mannheim, Rietz und Darmstadt‘ ziemlich 
gleichmäßig berücksichtigt sind. 

Dresden, Willers. 


13* 


Fr. A. Willers (o. Prof. a. d. Techn. Hochschule 
Dresden, Mathematische Maschinen 
und Instrumente. XII+ 318 S. mit 258 Abb. 
Berlin 1951. Akademie-Verlag. Preis geb. 34,— DM. 

Bei der lebhaften Entwicklung in der Konstruk- 
tion mathematischer Maschinen und Instrumente und 
bei dem regen Interesse, das insbesondere die großen 
Rechenautomaten und Integriergeräte mit ihren er- 
staunlichen Leistungen allenthalben finden, wird das 
Erscheinen einer zweiten Auflage des erstmalig 1943 
erschienenen Buches mit Freude begrüßt werden, 
Gegenüber der ersten Auflage ist der Inhalt beträcht- 
lich erweitert und dem gegenwärtigen Stand der Ent- 
wicklung angepaßt worden. Neu sind ein sehr klar ge- 
schriebener Abschnitt über Rechenautomaten, der 
eine gute Kenntnis aller wesentlichen Teile vermittelt. 
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Andere Erweiterungen betreffen Maschinen im Klein- 
format, Pantographen, Kurvenumwandler, Inver- 
soren, Schiebeplanimeter für Stieltjes-Integrale, Inte- 
griergetriebe als Differentiatoren, neue Analysatoren 
"und Integratoren sowie die mehrdimensionale Fourier- 
synthese. Schließlich trägt eine Umarbeitung des Ab- 
schnitts über Differentialgleichungsmaschinen den_ 
Fortschritten auf diesem Gebiet Rechnung, Der Um- 
fang des Literaturverzeichnisses ist von 489 auf 871 
Nummern angewachsen. vr; 5 . 

Das Werk und sein Verfasser sind so bekannt daß 
sich eine neue Empfehlung erübrigt. Erfreulich sind 
die gute Ausstattung, die instruktiven Zeichnungen 
und Bilder sowie ihre gute Reproduktion. 


Stuttgart. Günther Schulz. 


Die besprochenen und angezeigten Bilder sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.H. Ebert, Physikalisches Taschenbuch, 
VIIL-+ 5228. mit 147 Abb. Braunschweig 1951. 
Friedrich Vieweg & Sohn. Preis geb. 14,80. 


Dr. Ludwig Föppl (o. Prof. a. d. TH. München) und 
Dr. Ing. Gerhard Sonntag (Privatdoz. a. d. TH. Mün- 
chen), Tafeln und Tabellen zur Festigkeits- 
lehre. 226 S. mit 305 Abb. München 1951. Verlag 
R. Oldenbourg. Preis geb. 30,— DM. 


Dr. A. Pröll (Prof. i.R.a.d. Techn. Hochschule Han- 
nover, Grundlagen der Aeromechanik 
und Flugmechanik XVI+6128. mit 
278 Abb. Wien 1951. Springer-Verlag. Preis geb. 
48,— DM. E 


Prof, Dr. F. Henning, Von tiefen und 
hohen Temperaturen. Eine allgemein ver- 
ständliche Darstellung. 84 S. mit 44 Abb. Leipzig 
1951. B.G. Teubner Verlagsgesellschaft. Preis kart. 
3,40 DM. 


3. P. Den Hartog (Prof. am Massachusetts Institute 
of Teehnology und G.Mesmer (Prof. a. d. Techn. 
Hochschule Darmstadt), Mechanische Schwin- 
gungen. Zweite Auflage. XVI-+4278. mit 
299 Abb. Berlin-Göttingen-Heidelberg 1952. Springer- 
Verlag. Preis geb. 42,— DM. 


Th.Pöschl (Prof. a. d. Techn. Hochschule Karls- 
ruhe), Lehrbuch der technischen Me- 
chanikfür Ingenieure und Physiker. 
Zweiter Band:.Elementare Festigkeits- 
lehre. Zweite, umgearbeitete Auflage. VII+ 


je 


Internationale Mathematische Union. 
Vom 6. bis 8. März fand in Rom die erste General- 


- versammlung der International Mathematical Union 


statt. Dieser gehören z.Zt. folgende Länder an: 
Argentinien, Australien, Belgien, Canada, Cuba, Däne- 
mark, Deutschland, England, Finnland, Frankreich, 


Griechenland, Holland, Italien, Japan, Jugoslawien, . 


Norwe en,Österreich, Pakistan, Peru, Spanien, Schweiz, 
Vereinigte Staaten von Amerika. Außerdem nahmen 
an der Generalversammlung Beobachter aus Polen 
und Portugal teil.- Die deutschen Delegierten waren: 
Kamke-Tübingen, Knopp-Tübingen, H.L. Schmid- 
Berlin, Strubecker-Karlsruhe. Für die Zeit bis zum 
1.1.1955 besteht der geschäftsführende Ausschuß 
aus folgenden Personen: Stone-Chicago (Präsident), 


damm 19, Fernsprecher Sammelnummer; 4255 71. Postscheckkonto; Berlin 35 021. Bestell- dieses He x | 
1009/32/6. Die Zeitschrift für angewandte Mathematik und Mechanik erscheint monatlich, ne ee vierteljährlich. DM 15,— 
zuzüglich Bestellgeld, Einzelheft DM 6,—. Abbestellungen können nur bis 4 Wochen vor Quartalsende anerkannt werden, 5 j 
Demokratischen Republik, Satz und Druck: Druckerei Thomas Müntzer, Langensalza (K57 275 4022) (1). in German 


wird das folgende Quartalnoch geliefert. Veröffentlicht unter Lizenznummer 1207 des Amtes für Literatur und ara tree der Deut Yei, 
T FE 


nn 


2448. mit 159 Abb. Berlin-Göttingen-Heidelberg 
1952. Sprinßer Kur ae Preis geh. 16,50 DM, geb. 
19,50 DM. . 


Dr.-Ing. M. Strscheletzky, Hydrodynami- 
scheGrundlagenzurBerechnungder 
Schiffsschrauben. Mit einem Geleitwort 
von Prof. Dr. Betz. 257 8. mit 57 Abb. Karlsruhe 
1950. Verlag G. Braun. \ ‚ 


Prof. Dr. techn. Josef L. Krames (Mitglied der 
österr. Akademie der Wissenschaft in Wien) Dar- 
stellende und kinematische Geo- 
metrie für Maschinenbauer. Zweite, er- 
weiterte Auflage. VIII + 2678. mit 306 Abb. und 
2 Tafeln. Wien 1952. Franz Deuticke. Preis brosch. 
12,— DM; geb. 15,— DM. rs 


MAN Forschungsheft 1951. 95 S. m. 64 Abb. Augs- 
burg 1951. Verlag E. Kieser KG. 


Prof.L.Balser, EinführungindieKarten- 
lehre (Kartennetze). (Mathematisch-physikalische 
Bibliothek Reihe l Nr.81.) Zweite Aufl. 648. 
mit 50 Abb. Leipzig 1951. B.G. Teubner Verlags- 
gesellschaft. Preis 3,— DM. 


Dr.phil. F.A. Fischer (wissensch. Mitarbeiter im 


Fernmeldetechn. Zentralamt, Darmstadt). Grund- 
züge der Elektroakustik 
Bücherei) 162 S. mit 102 Abb. und 5 Tabellen. Ber- 


‚lin 1950. Fachverlag Schiele & Schön. Preis geb. 


14,— DM. j 


NACHRICHTEN  —- 


Borel-Paris (1. Vizepräsident), 
\ Vizepräsident), i 
yanaga-Tokio, Hodge-Cambridge ( and), Jessen- 
Kopenhagen. Für ern a un wurden 
Ausschüsse eingesetzt, darunter auch Internatio- 


Kamke-Tü 


nale Mathematische Unterrichtskommission (Imuk), 


(deutsches Mitglied Behnke-Münster). 
Tübingen 


Graz, Anlaßlich"einer Trinneronesfeire an rs 


den österr. Technischen Hochschulen vor 50 Jahren 
verliehene Promotionsrecht hat die Technische Hoch- 


schule Wien dem Prof. Dr. Karl Federhöfer, 
€ och 


rer: 


(Frequenz- 


Bompisni-Rom (Sekretär), 


E. Kamke, 2 > 


ar 
# 


